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We present the minimal supersymmetric standard model with general broken R-parity, focusing 
on minimal supergravity (mSUGRA). We discuss the origins of lepton number violation in super- 
symmetry. We have computed the full set of coupled one-loop renormalization group equations for 
the gauge couplings, the superpotential parameters and for all the soft supersymmetry breaking 
parameters. We provide analytic formulae for the scalar potential minimization conditions which 
may be iterated to arbitrary precision. We compute the low-energy spectrum of the superparticles 
and the neutrinos as a function of the small set of parameters at the unification scale in the general 
basis. Specializing to mSUGRA, we use the neutrino masses to set new bounds on the R-parity 
violating couplings. These bounds are up-to five orders of magnitude stricter than the previously 
existing ones. In addition, new bounds on the R-parity violating couplings are also derived demand- 
ing a non-tachyonic sneutrino spectrum. We investigate the nature of the lightest supersymmetric 
particle and find extensive regions in parameter space, where it is not the neutralino. This leads to 
a novel set of supersymmetric signatures, which we classify. 



I. INTRODUCTION 

The most widely studied supersymmetric scenario is 
the minimal supersymmetric standard model (MSSM) 
with conserved R-parity P, 0) 01 ■ "^^^ unification of the 
three Standard Model gauge couplings, gi, at the scale 
Mx = ©(lO^*^ GeV) iJI, is a strong indication that su- 
persymmetry (SUSY) is embedded in a unified model. In 
the simplest such model 111], SUSY breaking occurs in a 
hidden sector (decoupled from the Standard Model gauge 
interactions), and is communicated to our visible sector 
via gravity [3] ■ The scale of SUSY breaking in the visible 
sector is thus the Planck scale, Mp = 10^^ GeV. 

The large number of parameters in the MSSM is re- 
stricted by making well-motivated simplifying assump- 
tions at the unification scale. In the special case of the 
minimal supergravity model (mSUGRA), there are five 
parameters beyond those of the Standard Model: 



^0, Ml/2, Aq, tan/3, sgn{n) . 



(1) 



These are the universal scalar mass, Mq, gaugino mass, 
Ml/2, and trilinear scalar coupling, Aq, respectively, as 
well as the ratio of the Higgs vacuum expectation values 
(vev's), tan/3, and the sign of the bi- linear Higgs mixing 
parameter, /x. Given these 5 parameters at the unifi- 
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cation scale, we can predict the full mass spectrum as 
well as the couplings of the particles at the weak scale 
via the supersymmetric renormalization group equations 
(RGEs). This is the most widely used model for exten- 
sive phenomenological and experimental tests of super- 
symmetry. It is the purpose of this paper to create an 
analogous model in the case of supersymmetry with bro- 
ken R-parity (^p): the R-parity violating minimal super- 
gravity model (^p-mSUGRA). 

The mSUGRA model with universal boundary con- 
ditions was first extended to include bi-linear by 
Hempfling [^, focusing on the neutrino sector. A fur- 
ther detailed analysis in this framework was performed 
by Hirsch et al. Q. de Carlos and White were the first 
to go beyond bi-linear and consider the full set of 
-^p-couplings . However, they restricted themselves 
to the third generation Higgs- Yukawa couplings and used 
an approximate method to minimize the scalar potential. 
We detail below how we go beyond this work. 

We shall consider the chiral superfield particle content 



TT^ T °- TP TJ"^ TJ'^ 
n l-'l 1 t-^i 7 l^i, l^i, Hi, II2 ■ 



(2) 



Here i = 1,2,3 is a generation index, x = 1,2,3 and 
a = 1, 2 are SU{3) and SU{2) gauge indices, respectively. 
In supersymmetry, the lepton doublet superfields L° and 
the Higgs doublet superfield coupling to the down-like 
quarks. Hi, have the same gauge and Lorentz quantum 
numbers (This is an essential feature in our discussion 
below.). When appropriate we shall combine them into 
the chiral superfields £0=0,. ..,3 = ^^^1,2,3)- The 

gauge quantum numbers of the chiral superfields and the 
vector superfields are given in Table 
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TABLE I: The particle content of the mSUGRA _^p-model in terms of superfields and their decomposition into components 
with their SU{3)c x SU{2)l x U{1)y quantum numbers. x,X are SU(3) representation and generator indices, a, A are SU(2) 
representation and generator indices, a = 0, ..,3 is the family index of the lepton superfield, and i = 1, ..3 the usual family 
index of quarks, leptons and their superpartners. The fermionic components of the superfields are two component Weyl spinors. 



Chiral Superfields 



SU{3)c X SUi2)L X 



Components 



Df 

U't 

m 



(3,2,1) 
(3,1,1) 

(3,1,-i) 

(1,1,1) 
(1,2,1) 



UL 
dL 



UL 
dL 



d*R , dn 



h- 



e-Li 



eLa 



Sr , Sr 



h+ 
hi 



h+ 
hi 



h\ 



Vi 

eLi 



Vector Superfields 



SU{3)c X SU{2)l X t/(l), 



Components 



Vi 



(1,1,0) 
(1,3,0) 
(8,1,0) 



A. R-Parity Violation 

R-parity is defined as the discrete multiplicative sym- 
metry 



Rp - (-1 



(3) 



where S is the spin, B the baryon number and L the lep- 
ton number of the particle. All Standard Model particles 
including the two scalar Higgs doublets have Rp = +1, 
their superpartners have Rp = — 1 . When allowing for R- 
parity violation, the full renormalizable superpotential is 
given by 



(4) 



Here Ye,d,u are 3x3 matrices of Yukawa couphngs; 
Ayfc, X'ijj^, X'ljj^ are Yukawa couplings and Ki are mass- 
dimension one parameters, tab and exyz are the totally 
anti-symmetric tensors, with ei2 = 6123 = -1-1. The terms 
proportional to A, A', A", and ki violate R-parity explic- 
itly and it is their effect that we investigate in detail in 
this paper. The terms proportional to A" violate baryon- 
number, whereas the terms proportional to A, A', and Kj 
violate lepton number. Baryon- and lepton-number vio- 
lation can not be simultaneously present in the theory, 
otherwise the proton will decay rapidly [l^ ITsI . We dis- 
cuss in detail in Sect. ITU how this can be guaranteed. 

When extending mSUGRA to allow for R-parity viola- 
tion, the particle content remains the same but we have 
additional interactions in the superpotential, Eq. as 
well as the soft-breaking scalar potential (c./. Eq. H3U|l ). 



Thus within the ^p-mSUGRA the RGBs must be mod- 
ified. The running of the gauge couplings is only af- 
fected at the two-loop level and the effects have been 
discussed in Ref. Ref. |Q| also contains the 

two-loop RGEs for the superpotential parameters. Here 
we restrict ourselves to the one-loop RGEs. In order to 
fix the notation, we present the RGEs for the superpo- 
tential couplings as well as the gauge couplings in Ap- 
pendix ^ Due to the flavour indices the RGEs for the 
soft supersymmetry breaking terms are highly coupled 
to each other. In Appendix we discuss a very ele- 
gant method developed by Jack and Jones [l^ to derive 
the full set of RGEs for the soft-supersymmetry break- 
ing terms and apply it to the case of the ^p-mSUGRA. 
As we discuss, Jack and Jones' method is more easily 
implemented in a numerical computation. We also in- 
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dependently calculate the /3-functions of the theory by 
using the formulae from Ref. The resulting RGEs 

for the soft-supersymmetry breaking terms are given ex- 
plicitly in Appendix[ni We have checked that our results 
for the /3-functions in Appendices C and D are in full 
agreement. Furthermore, where relevant, they agree with 
previous (subsets of) results which have been computed 
by the standard method 8, 17]. 

Given the RGEs, we can compute the full model at 
the weak scale, including the mass spectrum and the 
couplings of all the particles as a function of our uni- 
fied scale (Mx) boundary conditions. In our numerical 
results for the ^p-mSUGRA, we extend the parameters 
given in Eq. by only one ^p-coupling. We thus have 

{A,A',A"}i, Mo, Ml/2, Ao, tan/3, sgn{fi) , (5) 

as our free parameters at AIx- {A, A', A"}i indicates that 
only one |?p-coupling is non-zero at Mx- We note that 
through the coupled RGEs many couplings can be non- 
zero at Mz and this is taken into account in the numerical 
implementation of our RGEs. 

Due to existing experimental bounds on the (A, A', A") 
[TsL H^ . the couplings are typically small and we thus 
expect the deviations from mSUGRA due to to be 
small. However, besides the RGEs discussed above there 
are four important aspects where there are significant 
changes and which we dicuss in detail in this paper: 
(i) the origin of lepton number violation, (ii) minimizing 
the scalar potential, (Hi) neutrino masses, and (iv) the 
nature of the lightest supersymmetric particle (LSP). 

(i) Since the discovery of neutrino oscillations, we 
know that lepton flavour is violated. If the ob- 
served neutrinos have Majorana masses, then lep- 
ton number is violated as well. In the ^p-MSSM, 
lepton number is naturally violated in the super- 
potential by the Yukawa couplings (A, A') as well 
as the mass terms Ki. In Sect. we discuss the 
origin of these terms in high energy unified theories 
and argue that they are just as well motivated as 
in the R-parity conserving case. For this we reana- 
lyze the seminal work on Z2 and Z^j discrete gauge 
symmetries by Ibanez and Ross j23|. We find a 
slightly different set of allowed operators, but the 
conclusions remain the same. 

We argue that within supergravity, with gravity 
mediated supersymmetry breaking, it is natural to 
have both = and = at the unification 
scale, Mx- This has not been taken into account 
in previous ^p-RGE studies. (Here Di is the to 
Ki corresponding soft supersymmetry breaking bi- 
linear term, c.f. Eq. I|3U|) .) This reduces the num- 
ber of parameters we must consider to the set given 
in Eq. 0. At the weak-scale however, in general 
Ki, Di 0, but these are then derived quantities. 

(ii) Since the lepton doublet superfields have the 
same gauge and Lorentz quantum numbers as the 



down- like Higgs doublet Hi, we effectively have a 
five Higgs doublet model for which we must mini- 
mize the scalar potential. Within our RGE frame- 
work, this must be done in a consistent approach 
while maintaining the value of tan /3 given at the 
weak scale and also obtaining the correct radiative 
electroweak symmetry breaking [23|. In Ref. 
(bi-linear R-parity violation) , points were tested to 
see if they minimize the potential for the case (in 
their notation) B — Di = Ag — 1. We have directly 
minimized the potential and do not make the lat- 
ter additional assumption. Instead, we determine 
B, Di via electroweak radiative breaking. If we ob- 
tain a point with radiative breaking of colour or 
electric charge we disregard it. We also go beypnd 
the numerical approximations made in Ref. to 
obtain the full result. The technical details of the 
iterative procedure are given in Sect. Hvl 

(Hi) Due to the coupled ^p-RGEs, a non-zero A or A' to- 
gether with /x(Mx)^0 will generate non-zero k^'s 
at the weak scale jMlSHSl- The k^'s lead to mixing 
between the neutrinos and neutralinos resulting in 
one non-zero neutrino mass at tree- level 24 .25...26.j . 
Thus one or more non- vanishing (A, A') at Mx will 
result in one massive neutrino at the weak-scale via 
the RGEs and the Ki. Requiring this neutrino to 
be less than the cosmological bound on the sum of 
the neutrino masses determined by the WMAP col- 
laboration 27] using their data combined with the 
2dFGRS data |3 

^m^,, < 0.71 eV, (6) 

thus gives a bound on the (A, A') at AIx- These 
bounds are determined in Sect. IVll Al and are very 
strict for the specific mSUGRA point SPSla, but 
are fairly sensitive to the precise choice of param- 
eters at Mx- The bounds are summarized in Ta- 
ble In Refs. 13, 113, it was argued that such 
bounds exist, however no explicit bounds were de- 
termined and the full flavour effects were also not 
considered. Here we present for the first time a 
complete analysis of the corresponding bounds. In 
a future publication we will address the possibil- 
ity of solving the atmospheric and solar neutrino 
problems within our framework. 

(iv) In the MSSM and mSUGRA the LSP is stable due 
to conserved R-parity. It can thus have a signifi- 
cant cosmological relic density Is^, |^ • Obser- 
vational bounds require the LSP to be charge and 
colour neutral 'si^ with a strong preference for the 
lightest neutralino, Xi- In -^p, LSP is not sta- 
ble and thus not constrained by the observational 
bounds on relic particles [32j ] . Therefore any super- 
symmetric particle can be the LSP: 

G, Xi, xf, ft=i,...,6, ^f=i,...fi, 1^3=1.2,3 , (7) 
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where x°, xt denote the Hghtest neutrahno and 
chargino, and qi, if, i>j, denote the right- and left- 
handed squarks, and charged sleptons as well as the 
left-handed sneutrinos, respectively. 

Depending on the nature of the LSP, the collider 
phenomenology will be completely different 34] . 
It is not feasible to study the full range of signa- 
tures resulting from the different possible LSPs in 
Eq. ^ or the different possible mass spectra. It 
is thus mandatory to have a well motivated mass 
spectrum, including the LSP, as in the MSSM and 
mSUGRA. Below in Sect. IVIII we determine the 
nature of the LSP as well as the rest of the mass 
spectrum as a function of our input parameters. In 
the no-scale supergravity models, we find signifi- 
cant ranges where the f is the LSP. In Sect. IVIIII 
we discuss the phenomenology of a f LSP. 

The case of a stau LSP has to our knowledge first 
been discussed in Ref. [s^, in the framework of 
third generation bi- linear R-parity violation. In 
Ref. |3a| the case of tri-linear was considered, 
with the focus on the comparison between charged 
Higgs and stau-LSP phenomenology. We go be- 
yond this to present a systematic analysis of all 
possible stau decays depending on the dominant 
|?p-coupling and classify the resulting signatures. 
For a recent analysis on charged slepton LSP decays 
in the presence of trilinear or bilinear ^p-couplings, 
see Ref. [S^. There only two-body decays are con- 
sidered and the parameters are restricted to the 
simultaneous solution of the solar and atmospheric 
neutrino problems. In Sect. IVIIII we present the 
general analysis. 

Very recently, in Ref. |32|| , the nature of the LSP 
in correlation with the neutrino properties is stud- 
ied in bi- linear ftp, i.e. the tri-linear couplings 
Kjk, ^ijk: Kjk to by hand. Since 

also the dependence on the supersymmetry break- 
ing parameters is not the focus of investigation, this 
work is complementary to our's. 

A stau LSP with R-parity conservation on the lab 
scale, i.e. the stau is stable in collider experiments, 
has been discussed in Ref. [s^- For completeness 
we mention that within R-parity conservation sev- 
eral authors have considered the case of a gluino 
LSP 113. 

B. Outline 

In Sect, present the motivation for supersym- 

metry with broken R-parity and discuss the possible ori- 
gins of baryon- and lepton-number violation. We focus in 
particular on the origin of the LiH2 mixing. In Sect. lIIll 
we present the full set of parameters and interactions in 
the mSUGRA model with broken R-parity, including the 
SUSY breaking parameters. In Sect. II Vl we discuss the 



radiative electroweak symmetry breaking including the 
full minimization of the Higgs potential. In Sect. [V] we 
determine the complete mass spectrum as a function of 
our parameters. In Sect. IVII we discuss the boundary 
conditions we impose at Mx and their numerical effect. 
In Sect. IV 111 we present our main results including the 
bounds we obtain on the _^p-Yukawa couplings from the 
WMAP constraint on the neutrino masses. In Sect. lVliTl 
we discuss the phenomenology of the stau LSP, classify- 
ing possible final state signatures at colliders and com- 
puting the stau decay length. We offer our summary and 
conclusions in Sect. IIXI 

We present two methods for computing these equations 
in the Appendices lAIBI and lO We present the complete 
set of RGEs at one- loop in Appendix [O In Appendix D 
we compute the four-body decay of the stau. 



II. ORIGINS OF LEPTON- AND 
BARYON-NUMBER VIOLATION 

In this section we investigate general aspects of the 
origin of baryon- and lepton-number violation in super- 
symmetry and thus the motivation for R-parity violation 
We then discuss in more detail the origin of the 
KiLiH2 terms in the context of only lepton-number vio- 
lation. In particular, for the following, we would like to 
know under what conditions after supersymmetry break- 
ing we can rotate away both the KiLiH2 terms and the 
corresponding soft breaking terms DiLiH2. 

A. Discrete Symmetries 

In the MSSM in terms of the resulting superpotential, 
R-parity is equivalent to requiring invariance under the 
discrete symmetry matter parity |12']. If instead we re- 
quire invariance under baryon-parity 

{Q,U,D) ^ -{Q,U,D), 
{L,E,Hi,H2) — > {L,E,Hi,H2) , 

we allow for the terms LLE, LQD, and LH2 in the su- 
perpotential, while maintaining a stable proton. Simi- 
larly lepton parity only allows for the U DD terms. Thus 
when allowing for a subset of R-parity violating interac- 
tions which ensure proton stability, we must employ a 
discrete symmetry which treats quark and lepton super- 
fields differently. In grand unified theories (GUTs) this 
is unnatural, as we discuss below. 

In string theories, we need not have a simple GUT 
gauge group. Thus models exist for both lepton- and 
baryon- number violation ,,41.] , and there is no preference 
for i?p-conservation or ^p. However, discrete symme- 
tries can be problematic when gravity is included. Un- 
less it is a remnant of a broken gauge symmetry the 
discrete symmetry will be broken by quantum gravity 
effects l43|. The requirement that the original gauge 
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symmetry be anomaly free, can be translated into a set 
of conditions on the charges of the discrete symmetry 
[2(tI . Considering the complete set of Z2 and Z3 dis- 
crete symmetries, and the particle content given in Ta- 
ble ^ only the Z2 symmetry R-parity, R2 , and the Z3 
symmetry B3 = R3L3 [4^ are discrete gauge anomaly- 
free '2^. i?3 is baryon-parity and allows for the interac- 
tions LLE, LQD and LH2 but prohibits UDD. 

This however does not completely solve the problem of 
proton decay. In supersymmetry there are also dangerous 
dimension-five operators which violate lepton- or baryon- 
number. The complete list is 





= [QQQL]f, 


O2 


= [UUDE]f, 


03 


= [QQQHi]f, 


O4 


= [QUEHi]f, 


05 


~ [LLH2H2]f , 


Oe 


= [LHiH2H2]f 


07 




Os 


= [HiHiE]F) , 


09 


= [QUL*]u, 


do 


= [QQD*]d, 



where we have dropped gauge and generation indices. 
The subscripts F, D refer to taking the F- or the D- 
term of the given product of superfields. We differ 
from Ref. 20] in that we have dropped the operator 
[H2H2e*]F), which vanishes identically and included the 
operator [QQd*]F). As in Ref. we have systemati- 
cally studied which Z2 or Z3 symmetry allows for which 
dangerous dimension-five operators. Our results are sum- 
marized in Table ^ We find some slight discrepancies 
with Ref. [231 . Furthermore, we have added the bilin- 
ear superpotential term KiLiH2 (K-term) not presented 
in Ref. |2(l |. As expected, the /i-term and the K-term go 
hand in hand in generalized baryon parity models (GBP) 
but the opposite is true for the generalized matter (GMP) 
or lepton (GLP) parity models: since the yit-term should, 
phenomenologically, be a non-zero parameter the GMP 
or GLP models containing the K-term are experimentally 
excluded. The requirement of neutrino masses excludes 
also the GMP and GLP models which do not allow for 
the AL = 2-term: LLH2H2. These models do not have 
any other source, within perturbation theory to incorpo- 
rate neutrino masses. From the models left, i.e. [GMP 
: i?2, GLP : L2, L3, GBP : R2L2, R3L3] only two can 
be induced from broken and anomaly free gauge symme- 
tries: these are the GMP: R2 (the usual R-parity case) 
and the GBP: B3 = i?3L3. 

Thus what we see from Table Ull is that although the 
MSSM R-parity is capable of eliminating the dimension- 
4 operators it is not capable of eliminating those of 
dimension- 5. Both dimension-4 and dimension- 5 baryon 
number violating operators are not allowed if the Z3- 
discrete symmetry B3 = R3L3 is imposed instead of the 
R-parity {R2) symmetry. In this article we study the 
phenomenology of the model based on the discrete Z3 
symmetry B3 = R3L3 |45|. 



B. Grand Unified Models 

In GUTs, quarks and leptons are in common multiplets 
and this simple approach does not suffice. We consider 
the case of the gauge groups SU{5) and 5*0(10) sepa- 
rately. 

1. SU(5) 

In SU{5) models, the trilinear and bi-linear R-parity 
violating terms are respectively given by 

h,jk *,*j^fc , k *«*5 , (10) 

where 5* is the 5* representation containing the D and 
L superfields, X is the 10 representation containing the 
Q, U, and E superfields and $5 is the Higgs superfield 
in the 5 representation, hijk are Yukawa couplings and 
ki dimension-one couplings. i,j,k are generation indices. 

Unless hijk ~ 10"^'^, this leads to unacceptably rapid 
proton decay. Thus this term must be forbidden by an 
additional symmetry. The generalization of matter parity 
where now and X change sign prohibits both terms in 
Eq. H10|) and guarantees that _^j,-terms arc not generated 
once SU{5) is broken. 

Alternative (discrete) symmetries can also be consid- 
ered. In Ref. 01, a discrete five-fold R-symmetry is con- 
structed which prohibits the terms in Eq. (|10|l . However, 
after breaking SU{5) SU{3) x 5f7(2) x C/(l) and in- 
tegrating out the heavy fields the operators ki vE'i'&s are 
generated, resulting in bi-linear R-parity violation. The 
size of the coupling depends on the vacuum expectation 
values of the large dimensional Higgs field representations 
which break SU(5). Similar symmetries can also be con- 
structed to obtain tri-linear R-parity violation. This was 
done in the case of "ffipped" SU{5) x U{1) in Ref. ^ and 
is easily transferred to the case of SU (5). The question of 
whether it is possible to obtain in GUTs with a large 
AL/AB hierarchy was also addressed in Ref. ^3 em- 
ploying a modified version of the minimal SU{5), where 
a built in Peccei-Quinn symmetry is broken at an inter- 
mediate scale. 



2. SO{W) 

In S'O(IO) GUTs B-L is a gauge symmetry and 
thus R-parity is conserved. Explicitly, the matter fields of 
a family are combined in a (spinorial) 16 representation 
and the operators 

R,jfc = 16, • 16j • 16fc , (11) 

are not 50(10) invariant. (Again, i,j,k are generation 
indices.) As in the SU{5) case, one would now expect to 
generate R-parity violating terms after breaking 50(10) 
and B-L. However, as shown in Ref. |49i] . surprisingly. 
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H1H2 


LH2 


QQQL 


UUDE 


QQQHi 


QUEHi 


LLH2H2 


LH1H2H2 


H2 HiE 


QUL* 


VD'E 


QQD* 


GMP: 
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R3AILI 
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TABLE II: In the left column we have the complete list of independent Z2 and Z3 discrete symmetries as in [2nl |. GMP, GLP, 
GBP denote generalized matter parity, lepton parity and baryon parity, respectively. In the top row we have the complete list 
of dimension-five operators which violate baryon- or lepton- number (c./. Eq. Q)- We have also included the operator H\H2- 
The symbol ^ denotes that the corresponding operator is allowed by that discrete symmetry. There are a few discrepancies 
compared to Ref. 2(1] . 



this strongly depends on the Higgs representations chosen 
to perform the breaking. 

If we include a 16_f/-Higgs representation to break 
50(10), as well as higher dimensional Higgs represen- 
tations we have the non-renormalizable operators 

Nyfc^f = 16, • 16j • IQk ■ 16h • G(H) , (12) 

where G(H) is a function of the higher dimensional Higgs 
representations. When the Higgs fields get vacuum ex- 
pectation values, 50(10) is broken and in general R- 
parity violating operators will be generated. The exact 
nature of the resulting R-parity violation depends on the 
employed Higgs fields and can be consistent with proton 
decay experiments 50] . 

Instead, we can explicitly exclude a IGh representa- 
tion and break 50(10) by a 126-Higgs representation 
[4^ . Since R^fe is an odd product of spinorial repre- 
sentations it is itself a spinorial representation. Without 
16ij there is now no spinorial Higgs representation and 
thus no 50(10) invariant combination 

Rz.fe • G'(H) , (13) 



where G'(H) is a general tensor product of Higgs repre- 
sentations. Thus after spontaneous symmetry breaking 
the operators Rijk can not be generated and there is no 
explicit R-parity violation in the theory. However, in 
principle R-parity can still be broken spontaneously with 
{v) 7^ or {v"^) 7^ 0, where v'^ is a right handed neutrino 
(which in this paper is only included in this discussion 
of 50(10)). With the absence of a 16// it was shown 
in Ref. [i^l that F-flatness at the GUT scale requires 
{v'^) = 0. This is also stable under the renormalization 
group equations. At the GUT scale we must also have 
{v) — 0, otherwise SU{2)l would be broken at Mqut- 
Similarly at the weak scale, we must demand {v) ~ 
in order to avoid an unobserved Majoran. Thus in this 
model R-parity is conserved at all energies and guaran- 
teed by a gauge symmetry |49j |. 



We conclude, that a priori there is no preference in 
supersymmetric GUTs for or against R-parity violation. 
Finally, we note in passing, that there exist few attempts 
in the literature to construct superstring models which 
accommodate the lepton number ^pCouplings j5l| . 
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C. Origin of the HiLiH2 Terms 

It is well known, that through a field redefinition of 
the Li and Hi fields, the terms in the superpotential 
Eq. Q can be rotated away at any scale |2J|. The full 
rotation matrix in the complex case was only given re- 
cently in Ref. |52l |. ^/ter supersymmetry breaking, how- 
ever, they can only be rotated away jointly with the cor- 
responding soft breaking terms DiLiH2, if and Di are 
aligned f2z,'2^. Even if they are aligned at a given scale, 
this alignment is not stable under the renornialization 
group equations |l,ll2ill3- However, if Ki and Di are 
aligned after supersymmetry breaking then we can choose 
a basis where /t^ = = at the supersymmetry break- 
ing scale. At the electroweak scale, we then have a pre- 
diction for both Ki and Di through the renormalization 
group equations (RGEs), given the initial choice of basis. 
We are thus interested in the conditions for alignment 
after supersymmetry breaking in various unification sce- 
narios, in order to predict Ki{Mz) and Di{Mz)- 

We first consider the general superpotential of Eq. (0}, 
restricted for the case /i = Ki = 0). It is invariant under 
a discrete R-symmetry j5^ . where the chiral superfields 
have the following R-quantum numbers [s^ . 











Di 


Hi 


H2 





-2 


-1 


-1 


-1 









The vector superfields have zero charge. Each term in 
the superpotential must have R-charge -2, which is can- 
celed by the charges of the Grassman coordinates. Thus 
all tri-linear terms except UDD are allowed. Note, that 
since this is an R-symmetry the fermionic components of 
the chiral and vector superfields have a different charge 
than the superfield. In particular, the R-parity even com- 
ponents of the chiral superfields have the quantum num- 
bers of conventional lepton-number. With this somewhat 
unusual symmetry we have ensured Icpton number con- 
servation for the SM fields [sdj. 

However, the phenomenology of this superpotential is 
unacceptable. Below we show that if fi, Ki, B, Di = 0, 
the CP-odd Higgs boson mass, rn^ = and the lightest 
chargino mass M-± - 0(30 Gey), both in disagreement 
with observation, — due to the Peccei-Quinn sym- 
metry of the superpotential. We thus demand Ki, /i 7^ 0, 
in order to get consistent SU{2) x U{1) breaking and 
a sufficiently heavy chargino. This in turn introduces 
lepton-number violation for the low-energy SM fields. 

The parameters Ki and fi are dimensionful and in prin- 
ciple present before supersymmetry breaking. The only 
mass scale in the theory is the Planck scale (Mp), and 
we thus expect Ki, fi — 0{Mp). Experiment requires 
fi = 0{Mz) and Ki <C Mz- (The latter strict require- 
ment is due to neutrino masses, as we discuss in detail 
below.) This is the well-known //-problem [s^ . modified 
by the presence of the Ki . In the following, we discuss the 
origin of the weak-scale fi and Ki terms and their corre- 
sponding soft terms. We can then determine under what 



conditions the Ki and Di can be simultaneously rotated 
away at the unification scale. We begin by discussing 
supergravity theories where there are several proposed 
solutions to the //-problem [s^ |5^ [s^ ■ We review 
them here in the light of the additional Ki terms. 



D. Supergravity 

We consider a set of real scalar fields Zi for the hid- 
den sector and a set ya for the observable sector Qj. 
Collectively we denote them Za- The supergravity La- 
grangian depends only on the dimcnsionless scalar func- 
tion (Kahler potential) |^ 

g{z,,z";ya,yn = -d{z„z";ya,yn/M], 

~\og{f{z,;ya)/Ml). (14) 

Here d determines the Kahler metric and / is the super- 
potential, which is a holomorphic function. The scalar 
potential is given by 

1 



V = -Af4exp(-^) [3 + gA(g-i):^g^] 

( d{ZA,Z^*) \ 

= [ mi ' " 



(15) 



X [{d-')iF^^FB--if\Z^*)f{ZA)/Ml] + \d^D^ 



Here = dG/dZA, and 



Fa = 



dZAdZB* 
_ df(ZA) , ^^_2dd{ZA,Z^*) 



OZa 



M- 



OZa 



(16) 

f{ZA). (17) 



and £)" is the auxiliary field of the vector superfield. The 
derivatives of are defined analogously. 

The most general form of the low-energy scalar poten- 
tial after supersymmetry breaking is jH3 

+m3/2[%) + /i^(2/)] + ^i5"i?" . (18) 

Here g{ya) is the superpotential for the low- energy fields 
derived from /{Za) and m3/2 is the gravitino mass. The 
first and the last terms are the usual F- and Z?-term con- 
tributions to the scalar potential. The second and third 
terms arise from supersymmetry breaking. The general 
constant matrix Sab has in principle arbitrary entries, i.e. 
the soft scalar masses can be non-universal. 

h{y) is a superpotential, i.e. a holomorphic function of 
the ya- In the renormalizable case, it is at most trilinear 
in the fields ya and contains the supersymmetry breaking 
A and i?-terms ^6^. g{y) and h{y) are superpotentials of 
the same fields and due to gauge invariance thus contain 
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the same terms. However in general, the coefBcients are 
independent and thus in particular the A- and B-terms 
need not be proportional to the corresponding terms in 
g{y)- But if the superpotential satisfies 



(19) 



then the soft-breaking term h{y) is a linear combination 
of the superpotential gijj) and yadg{ya) / dya 62] and 
thus each term is proportional to the corresponding term 
in g{y). The condition (|19l) is quite natural. If the Zi 
all transform non-trivially under only the hidden-sector 
gauge group and the ya transform non-trivially only un- 
der the observable sector gauge group, then combined 
with the requirement of renormalizability we obtain the 
condition ((T^ . 

We now consider the observable sector superpotential 
given in Eq. If our superpotential at the unification 
scale satisfies Eq. ()19|l . the Di will be aligned with the 
Ki after supersymmetry breaking and they can be simul- 
taneously rotated away. Or looked at differently: before 
supersymmetry breaking we can always rotate the fields 
such that Ki — 0. If we then break supersymmetry 
at this scale, while obeying Eq. ^1 we automatically ob- 
tain Z?i = as well, since the coefficients in h{ya) are 
proportional to those in g{ya)- Thus in the case of a 
renormalizable superpotential we expect universal A and 
B terms and thus an alignment of and Di at the uni- 
fication scale. 



E. Implementing a Solution to the /^-Problem 

The most widely discussed solution to the /i-problem is 
to prohibit the /ii?ii?2 in the superpotential via a symme- 
try, for example an R-symmetry, and instead introduce 
a non-renormalizable term into the Kahler potential, 
which results in the /Lt-term after supersymmetry break- 
ing. By using the mass scale inherent in supersymme- 
try breaking one then obtains /x = 0{AIz)- This was 
first proposed by Kim and Nilles who introduced the 
non-renormalizable term into the superpotential /. The 
R-symmetry was global and the resulting axion w as p he- 
nomenologically acceptable. Giudice and Masiero [53 in- 
troduced a non-holomorphic term into the Kahler metric 
function d instead, also invoking an R-symmetry to pro- 
hibit terms in the superpotential. The details of the axion 
were not considered. In certain cases the two mechanisms 
are equivalent 53]. In the following, we briefly consider 
the implications of Ref. |58j for the Ki terms and extend 
this to Ref. ^]. 

In the context of R-parity violation we have both a /i 
and a Ki problem. As an example, we introduce the fol- 
lowing non-renormalizable terms into the superpotential, 

/' = 4r{aziZ2HiH2 + b^Z3ZiL,H2), (20) 
M p 

assuming them to be invariant under the symmetries of 
the model. In general, we could have higher powers of 



the Zi. If the Peccei-Quinn '64*] charges which prohibit 
the bilinear terms in the superpotential are lepton-flavour 
blind but distinguish Hi and Li then we would expect 
the general form shown above, a, bi are dimensionless 
constants. Due to the independent fields Zi we can not 
rotate away the bi terms. After supersymmetry breaking 
we get 



(^l)(^2) 

Mp 
Mp 



= 0{Mz) 
- 0{Mz) 



(21) 
(22) 



If the fields Zi are hidden-sector fields and /' mixes the 
hidden and observable sectors then the soft supersymme- 
try bilinears are in general not aligned with the Ki since 
there are now the additional terms 



1 

Mp 



dfi 
dz4 



dh 
dz3 



(23) 



which have independent coefficients from the purely hid- 
den sector. Here we have made use of the hidden-sector 
function /i of Eq. H19|) . The resulting Ki terms are still 
0{Mz)- If 9/i/9zi=i.2,3.4 = then we have alignment. 

Alternatively, the Peccei-Quinn charges can be such 
that Hi has the same charge as the Li. This is exactly the 
case of the Bt, — R3L3 discrete symmetry we discussed in 
some detail in section Hi Al and follow in this article. The 
charge of Hi and the Li under this symmetry is —2/3. 
In this case Z1Z2 — Z3Z4 in Eq. (|20() and the Ki terms 
can be rotated away before supersymmetry breaking. No 
Di soft terms are generated in supersymmetry breaking 
then and we have Ki — Di = at the high scale. 

We conclude that it is possible to have alignment of the 
bilinear terms at the supersymmetry breaking scale but 
not necessary. The eventual answer will depend on the 
underlying unified theory. We shall assume that we can 
rotate away the Ki terms before supersymmetry breaking. 



III. THE MINIMAL R-PARITY VIOLATING 
SUPERSYMMETRIC STANDARD MODEL 

The model we consider has the particle content given in 
Table m and the superpotential given in Eq. 10}. Within 
this superpotential we shall make the assumption that at 
the unification scale, Mx — 10^^ GeV, the terms KiLiH2 
have been rotated to zero. For real parameters the or- 
thogonal rotation on the fields which accomplishes 
this is given by 



and explicitly in components 



(Hi\ 




( 


-S3 





^ 






Li 




C2S3 


C2C3 


-S2 







^'1 


L2 




C1S2S3 


C1S2C3 


C1C2 


-si 




^'2 






\ S1S2S3 


S1S2C3 


S1C2 


Cl ) 







(24) 



,(25) 
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where Ci = cos 0i and Si = sin 9i , and 



ci = 



C2 = 



H2 



Si 



S2 = 



K3 



(26) 



C3 



S3 



Here we have introduced the notation — nf. The 
more general case of complex parameters is given in 
Ref. we shall restrict ourselves to real parameters 
here. After the above field redefinition, the only remain- 
ing superficld bi-linear term is 



(27) 



with /i" = ^y fj? + and H[ = C'q. This will be our 
starting bilinear superpotential term at Mx in our RGE 



studies below. 

The RGEs for the Ki are given by (see Appendix A) 



16^2_^. = ^^j^l + ^j^- + (28) 

where at one-loop the anomalous dimension mixing Li 
and Hi is given by 



ijTcl I E)jk 



(29) 



with a summation over j, k implied. (The remaining 
anomalous dimensions are given in Appendix 0) There- 
fore, given /i 7^ at Mx and a non-zero A or A',we will 
in general generate a non-zero Ki{Mz) H, 0, H^, l23| . 
Below we discuss special exceptional cases where this is 
not the case. 

In order to fix all the parameters we also need to know 
the general soft supersymmetry breaking Lagrangian 
which we denote 



2^xyzll'ijk^i ^j^k 



H.c 



(30) 



Here, F e [Q, U, D, E, £] denote the scalar component 
of the corresponding chiral superfield. rrip^ are the soft- 
breaking scalar masses. Note that these are 3x3 matrices 
for the squarks and for the lepton singlets. However, 
(m^^)^^ is a 4 X 4 matrix. (hc/)y, h^^k, h'^-^., and h'f^^. 

as well as ba — {B, Di) are the soft breaking trilinear and 
bilinear terms, respectively. 



The RGEs for the Di are given at one-loop by 



167r 



2^(71 



'Hi 



(7i) 



(7i)f>' 



(31) 



with the anomalous dimensions (7) and the functions 
(71) defined in Appendices IBl and O respectively. These 
RGEs are clearly distinct from those for Ki, above. It 
is thus clear that given Ki{Mx) = Di{Mx) = we 
will lose alignment between the two at the electroweak 
scale H,Cl4, 22, 23]. In order to describe the weak-scale 
physics, we thus require the full set of parameters given 
in Eqs. Q and (|2DJ. 



IV. ELECTROWEAK SYMMETRY BREAKING 



The full scalar potential is given by 

VscALAR, = VsuSY + VsoFT 



(32) 



with the supersymmetric F-term and D-term scalar po- 
tential given by [6^ 

VsusY = Vp + Vd = 



E 



dW 



3 2 

El 

1=1 



A 



, (33) 



respectively and VgoFT = — 'Csoft- In Eq. the 
fields ^m,n denote the scalar fields in the theory, 3^=1,2,3 
are the gauge couphngs with gi for [/(l)y, g2 for SU{2)l 
and ^3 for the SU (3)c gauge group. In order to simplify 

the expressions, we shall use the coupling g = \J\gi- 
TO, n.. and A, B... are representation and gauge generator 
indices, respectively. The explicit expressions for Vp and 
Vd can be found in Ref. 66] . 
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In the following, we shall focus on the complex neutral 
scalar fields: i^a = (^i, 1^^=1,2. 3)- For these the scalar 



potential is given by 



Neutral 



9l) 



AV, 



(34) 



r 



where A V denotes higher order corrections • In order 
to minimize this potential, it is convenient to write the 
complex neutral scalar fields in terms of CP-even, 0:2, Tq,, 
and CP-odd ?/2, real field fluctuations 



X2 



W2 : 



(35) 
(36) 



At the minimum the scalar fields thus take on the values 
{X2) = Vu, (ra) = Va, {va = {vd, vi , V2, V3)) and (j/2) = 



(ta) — 0. The minimization conditions for VNcutrai can 
be written as, 



dV 



Neutral 



dx2 



dV 



Neutral 



dra 



0, (37) 



where "min" refers to setting the scalar fields to their 
values at the minimum. We then derive the following 
five minimization conditions, where a, (3 = 0, 1, 2, 3 and 
there is an implied sum over repeated indices 



J 



3?e 



+52) \vu\-Kr + 



1 dAV 



2 dv(3 



= 



(38) 



Here SRe denotes the real value and we have written 
{dAV/dra)Un as dAV/dv^ and (aAV/5a;2)Uin as 
dAY/dvu- Next, we solve this system of equations. We 
start by defining 



vl = sin^/3(t,2~^t,2 



(42) 



and 



V = v: 



tan p = — , 

Vd 



^ 9 A f2 

l + v'd+T.-'- — 



2 ' 

52 



(39) 



(40) 



with Vi being the three sneutrino vev's. The advantage 
of using the definition given in Eqs. H39II40|) is that tan/3 
is the same in the R-parity conserved (RPC) and ^p- 
models. This facilitates the direct comparison, in partic- 
ular when Vi/v <C 1. 



where in our convention w = 174 GeV. Then the vev's Vd 
and Vu can be written 



Vj = COS' 



^ 1=1 ^ 



(41) 



Using these definitions and the notation (vf) = J^i 
the five minimization conditions in Eq. (|38|l can be writ- 
ten as (again there is an implied sum over repeated in- 
dices.) 



(wl/i + tJ-'^)vd + 



v,^ Bvu + l-Mlcos2Pvd + lig^ +gl)sm^ Pvd{vf) + l^^ = 0, (43) 
2 2 2 OVd 
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(m^^^J Vd+ (mL^)ji + K*K, t;j -At^« +iM|cos2/3w, +5|)sin2/3t;, = 0, (44) 

(to^^ + + \k^\'^)vu - Bvd - DiVi - ^Mf cos2(3vu ' \i9^ + fll) sin^ P Vu («■ ) + = • (45) 



r 



In order to solve the above equations, we first derive /i 
in terms of Wd, and Vi from Eqs. (|43|l and H45(l . It is 
obtained after solving the quadratic equation 



A^t^ + + r = 



(46) 



r = 



tan^/? 



+ -M|(tan2/3-l) . (48) 



with 



A EE tan^/J-l, B = — 

I'd 



The solution to Eq. (|46l) can be written in a more familiar 
(47) form, 



m2 



tan^ (3-1 



tan^/J 



(49) 



r 



We recover the familiar minimization condition 69] in 
the RPC limit k^, w^, A, (m| ^ ) ^ 0. 

Eq. H46|l or equivalently Eq. (|49|l . has two solutions for 
the parameter fj, : fi > and < 0. We thus retain 
the sign of /it as a free parameter. Furthermore, the fac- 
tor that multiphes the /i parameter in Eq. is 
small since, as we show below, Vi <C Vd to obtain a small 
neutrino mass, nii, ^ 0{eV). 

We can now express B in terms of /i, Vu, Vd, Vi from 
Eqs. 63 and ESJ 



B = 



sin 2/? 



2|Ai|' 



_ ~ 

Vd Vu 



(50) 



where in both Eqs. (|48l) and H5U|) we have introduced the 
simplifying notation 

1 dAV 



2 Vu dvu 
_2 _ 2 1 9AV 
^ 2vd dvd 
Eq. H44(l can now be cast in the form, 



(51) 
(52) 



Vd 



where 



(Ml), 



1 dAV 

^2^ 



(niL )ji + KiK* + -Mz cos 2/3 4 



(53) 



(5^ 



.91) 



sin^ /3 (w^ - vl 



(54) 



Here we outline the iterative numerical procedure we 
follow to obtain the minimum of the potential for a given 
value of tan/?. 

1. We start in the RPC limit with Vi = and thus 
obtain from Eqs. H41II42|I initial values for Vu and 
Vd (in terms of tan (3) . 

2. We solve Eqs. (gHJ (or and also first in 
the RPC limit, v, = 0,k, = D, = (m^^. ) = 0, 

and thus obtain initial values for /i and B. 

3. We treat w„, Vd and /x, B as known and solve the 
system of Eqs. (|53f) in terms of the Vi. This system 
is linear and a lengthy analytical expression of the 
solution exists. 

4. We return to the first step and compute the cor- 
rected values of Wu, Vd including the ViS using 
Eqs. l(lll|121). The reader should note that tan/3 = 
Vu/vd remains exactly the same as in the R-parity 
conserving MSSM case (see Eqs.l^ES- This is 
the advantage of this formulation-developed for the 
first time in Ref. [63 -and is used throughout this 
paper. In our calculation, we include the full one 
loop corrections and the dominant two loop ones 
as they have been calculated in the RPC case in 
Ref. Ipti but not R-parity violating loop correc- 
tions |67| . 

5. We repeat the second step but use the non-zero 
values of Vi as well as the newly computed values 
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of Vu, Vd- At this point we now also include the 
non-zero values of Kj, Dj. The latter could have 
been included from the beginning but it is compu- 
tationally more convenient to do this in the second 
iteration. 



where again v'^ 



The photon and the 



gluons are of course massless. The reader should note 
the participation of the sneutrino vev's Vi in the masses 
of the Z- and ly^-gauge bosons. 



6. We now iterate the procedure until convergence of 
Vu, Vd,Vi is reached. 

We have explicitly checked that our iteration procedure 
is very robust and for all the initial parameters we dis- 
play in our numerical results, we have found the iteration 
procedure to converge. 

Finally, it is well known that the MSSM provides 
a mechanism of breaking radiatively the electroweak 
SU{2)l X U{1)y symmetry down to C/(l)om i^J- Elec- 
troweak symmetry breaking in the MSSM occurs when 
mjj^ < in Eq. This is indeed realized in the MSSM 
since is driven to negative values by the large top 
Yukawa coupling once we employ the RGEs. As we see 
from Eq. (|(]18|l the -^p-couplings do not affect directly the 
"running" of . However, they do affect the running 
of rnj^^ in Eq. ljC17l) through the mixed wave function 
Hi — Li. These corrections turn out to be small, since 
'^'l Hi '^^ small, in the minimal supergravity scenario we 
assume in this article. Concluding, the radiative elec- 
troweak symmetry breaking in the ^p-case works in ex- 
actly the same way as in the RFC case. 



PARTICLE AND SUPERPARTICLE MASSES 



In the literature, it is common to make a specific ba- 
sis choice for the CP-even neutral scalar fields I'a, in 
particular the basis where only v^, Vd ^ and Vi — 0. 
We shall present our results for particle and superpar- 
ticle masses in the generic basis, where all vev's can be 
non-zero, w„, Vd, Vi ^0. We shall strictly follow the con- 
ventions of Grossman and Haber [b^ which in the R- 
parity conserved limit coincide with those of Haber and 
Kane 0. We list in turn the mass matrices and show 
how they depend on our basic parameters, as well as the 
minimum of the potential determined in the previous sec- 
tion. It is then straightforward for the reader to choose 
his/her favorite basis or to work with the basis indepen- 
dent spectrum given below. 



B. CP-Even Higgs-Sneutrino Masses 

From Eq. (|34|l . we see that after electroweak symmetry 
breaking, the sncutrinos, Vi, mix with the Higgs bosons 
hi = Dq. If CF is conserved, the mass eigenstates sep- 
arate into CF-even and CF-odd states. Following Gross- 
man and Haber 66], let us denote with the CF- 
even (CF-odd) sneutrino mass eigenstates. If R-parity is 
broken, the mass of !>+ is in general different from the 
mass of f'-, i.e. there is a sneutrino, anti-sneutrino mass 
splitting. The CF-even Higgs-sneutrino mass eigenstates 
are denoted by , H'^ , , where the mass M/iO < Alfjo . 
They are obtained in the generic basis after the diago- 
nalization of a 5 x 5 mass matrix 




where 



CP— even 



(57) 



(58) 



/ bqVa I (.9^ + 32) ,.2 isL+jl),, 

H 2 2 ^"^'5 



with 



(9^ + 92) f 2^ , (9^+5!) 



{9' +91) 



{vl - vl)5^p , (59) 



A. Gauge Boson Masses 



and where = ^ v'i. Recall that ba = (S, Di). 



C. CP-Odd Higgs- Antisneutrino Masses 

The CF-odd Higgs-sneutrino mass eigenstates A, 
(and the massless Goldstone boson in the unitary gauge) 
are obtained in the generic basis after the diagonalization 
of a 5 X 5 mass matrix 



For completeness and in order to fix our notation be- 
low, we write here the masses of the Z and gauge 
bosons. 



Ml 



2 / 2 I 2 \ 

92 («u +«a) : 

{9'+9l){vl + vl) 



(55) 
(56) 



where 



C 



'^(y2,2;7)-^cp-odd ( 
2 \ys 



CP-odd 



Vu 



(60) 



(61) 
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For one generation, we obtain two nonzero eigenvalues 
with the eigenstates identified as the sneutrino and the 
CP-odd Higgs, respectively, 



and A^cp-odd given above. Eq. H64|l leads to the follow- 
ing Higgs mass sum rule in the _^p-scenario. 




,(62) 



.(63) 



Here mj, is the one generation version of Eq. (|59|l . No- 
tice the tan (3 enhancement (reduction) of the sneutrino 
(Higgs) mass is due exclusively to an R-parity violat- 
ing contribution. For Di — > we have mj,_ = mj, and 
Mio = as it should be. 

The generalization of the Higgs mass sum rule M"^ 



in the RFC case is written here as: 



Ml 



Tr(Al2cp-o, 



•dd) 



(64) 



This is easily verified from the matrix forms oi 



Mto ^ Mho + 



Ml 



i=l 



(65) 



This sum rule is valid only at tree level and is altered 
by radiative corrections. If the heavy Higgs mass states 
and are degenerate and also the sneutrino anti- 
sneutrino mass difference is small then the light Higgs 
boson mass would be very close to the Z-boson mass. 



D. Charged Higgs Bosons-Sleptons 

The charged Higgs bosons mix with the charged slep- 
tons. 



C = -(/i2 :eL.,,eR J Charged 



/ h+ 



(66) 



CP— even 



In the basis independent notation, the 8x8 mass matrix 
is given by 



■^Charged 



I 



(67) 



with 



(m^)ii = mjj^ + I^Qp , 



{D)ik 



(68) 
(69) 

(70) 
(71) 

\{9l - 9^){vl ~ vl)^Sj + ^givjvs , (72) 
1 



(73) 



The remaining parameters are given in Eqs. Q and 
(|30|l . Upon diagonalization of the mass matrix (|67|l . we 
obtain the mass eigenstates : G^, i?''^, e,;=i^...^6- K is 
not hard to prove that the determinant of H67|) is zero 



and the Goldstone boson corresponds to the eigenvector 
(-w„,Wq, 0,0,0). 



E. S quarks 

1. Down Squarks 

The down squark mass eigenstates di,i — 1, . . . , 6 are 
given by diagonalizing the following mass matrix 



dR,+ 



(74) 



where in the {di-, d^j.^g} basis we have 
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ML 



(75) 



The * denotes the complex conjugate of the transposed 
matrix element, i.e. in the above case (■A^Down)!^- 



given in the {uL^,Ufi^^^} basis 



2. Up Squarks 



C = 



-(^L.'^fl.+3)-^&p 



(76) 



The up squark mass eigenstates Ui, i — 1, . . . , 6 are 
determined by diagonalizing the following mass matrix 



where 



9l 



(77) 



F. Quarks 

The down quark masses are given by, 

(i^dh = X'^.jVa , (78) 
and the up quark masses are 

(m„)ij = {Yu)ijVu , (79) 
and the coupling constants are defined in Eq. I^J. 

G. Neutrinos-Neutralinos 



nos. The neutrino-neutralino mass matrix (7x7 for three 
generations of neutrinos) in the {—iB, —iyv''^\ h^, i^a) ba- 
sis is given by 



C 



'-(-iB,-iW'-^\hliyc.)M^ 



-iB \ 

\ ^0 J 



, (80) 



The neutrinos mix with the neutralinos resulting in one 
massive neutrino at tree level and four massive neutrali- where ITC 
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Ml 



Mzs 



zsw 



-Mzsw 



Mzsw 



Mo 



-Mzcw 



Mzcw 



Mz cw Mz cw 



(81) 



with M^ given in Eq. H5t)|l and sw = sin dw , is the elec- with 
troweak mixing angle. The matrix H81|l has five non-zero 
eigenvalues, i.e. four neutralinos and one neutrino. We 
denote the mass eigenstates which are obtained upon di- 
agonalization of the matrix as: Xi 4j ^i=i....,3j with the 



masses M^o < M^o < M^o < M^o . 

^ 

Since Mi, M2, M^ S> Vi, the matrix Eq. l|HT|l is sugges- 
tive of the well known sea-saw formula. 







where M^ is the 4x4 neutralino mass matrix with mass 

eigenvalues typically M^. - 0(10 GeV) [nl- The off- 
diagonal entry m is a 3 x 4 matrix with entries of order 
gvi, or Ki. In Sect. IVlTI we show and below we estimate 
that Ki ~ 0(1 MeV) and thus m ^ M^. The analogy 
with the Majorana see-saw mechanism is then obvious 
under the replacements 



Mx = MsusY 



MMaj , 
Moirac 



(83) 



In addition, the 3x3 zero mass matrix in Eq. H82|l can be 
filled by finite, loop low energy threshold corrections in 
the ^p-MSSM as opposed to possible Higgs triplet con- 
tributions in other neutrino mass models. Therefore neu- 
trino masses will roughly be given by 



7 1 

an 



Ms 



USY 



Ms 



1 eV 



USY 



(84) 



For the last inequality, we have imposed the bound from 
WMAP in Eq. Bearing in mind possible accidental 
cancellations (see below) we obtain 



~lMeV for MgusY^lTeV. (85) 



A complete calculation of the one neutrino mass eigen- 
value at tree level reads HI, Hal 



(86) 



2vuVd{Mxgl -f M2g2) - 2nMiM2 ' 



Ai^Vi- Vd — 



(87) 



A redefinition of the phases of the gaugino fields B and 
W together with the gaugino universality assumption 
Ml — M2 = Ml/2: can make Mi and M2 real and posi- 
tive and so the numerator of Eq. (|86|l cannot be fine tuned 
to zero (provided /i > 0(100 GeV)). According to the 
universality assumption, the 1-loop unification gaugino 



(82) 

masses at the electroweak scale are, Mi 



3 ai 



-M 



^1/2 



and M2 = Tx^Mi/2, where acuT = gGUT/47r ~ 0.041 

"gut 



is the grand unified coupling constant. Taking into ac- 
count that VuVd ^ mMi/2, which we find in our numerical 
results below, we arrive with an excellent approximation 
at a simple formula for the tree level neutrino mass: 

_ 167raGUT ELi^i 



5 Ml/2 

This implies - 1 MeV for M1/2 - 1 TeV. One can 
obtain a small even with Vi ^ Ki ^ v but that requires 
a cancellation of 1 part in 10^. So the question arises 
how one can naturally obtain a small A^, i.e. Vi ^ Ki ^ 
0(1 MeV)? We will come to this point in Sect.lVlIl 



H. Leptons-Charginos 

The charged leptons mix with the charginos. The La- 
grangian contains the (5x5 for three generations of lep- 
tons) lepton-chargino mass matrix as ^6 

I 



C 



\ 




(89) 



where the mass eigenstates Xi2'^ = is.,fi,T) are given 
upon the diagonalization of the matrix 



Ml 



M2 



g2Vu 



Ok 



g2Va fJ'a ^(iakVp 



(90) 
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VI. BOUNDARY CONDITIONS AT Mx 

Due to the large number of parameters in the super- 
symmetry breaking sector (c./. Eq. (|3(J|l ). we shaU focus 
on the case of minimal supergravity models. These have 
a much simplified structure at the high scale, which we 
assume here to be the unification scale of the gauge cou- 
plings, Mx = Mgut = 0(10^*^). At this scale, the soft 
SUSY breaking scalar masses have a common value, Mq: 

m^iMx) = ms(Mx) = mj(Mx) = 
miiMx) = mg(Afx) = Moi (91) 

mnAMx) = mH,{Mx) = Mo, (92) 

where 1 is the 3x3 unit matrix in flavour space. Mo- 
tivated by the discussion of Sect. Ill, we shall assume 
that we can rotate away the Ki terms before supersym- 
metry breaking and no Di or (m? ^ ) terms are gener- 
ated through supersymmetry breaking at the unification 
scale Mx, 

K,{Mx)^0, D,{Mx) = {ml^^JiMx) = . (93) 

At the scale Mx, we shall assume one non-zero ^p- 
coupling at a time, i.e. one coupling from: 

KMMx)^0 , A^fe(Afx)y^O , X';^,{Mx)^0.{M) 

Due to the CKM quark mixing, the A' RGEs are coupled. 
Thus in the case of a single A'(Mx) 7^ we will have 
more than one X'{Mz) ^ at the weak scale. mSUGRA 
assumptions lead to the same prefactors, Aq of the su- 
persymmetry breaking trilinear couplings hijk = AoYijk'. 

Au(Mx) = Au{Mx) - Ae{Mx) - 

Aa(Mx) = Ay{Mx)=Ax„{Mx) = Aoi. (95) 

A common mass, Afi/2 for the gauginos completes the 
mSUGRA boundary conditions at Mx, 

Ml (Mx ) = M2 {Mx ) = M3 {Mx) = Ml/2 ■ (96) 

No assumption for quark or lepton Yukawa unification 
has been made in our analysis. We thus have the six 
parameters : 

Ao, Mo, Ml/2, tan/3, sgn{^i), {A,A',A"}i. (97) 

When determining the mass spectrum, in order to further 
simplify the number of input parameters we will restrict 
ourselves to a particular supergravity scenario called "no- 
scale" supergravity [t^I . This scenario predicts a definite 
relation between Aq and Mq namely 

Aq = Mo = GeV . (98) 

The "no-scale" scenario, the simplest mSUGRA sce- 
nario, is experimentally excluded in the RFC case, but 



as we show below allowed in the ^p-case. Our re- 
sults for the bounds on the .^p-couplings from neutrino 
masses should be unaffected by this assumption provided 
{Mq, \Ao\)/Mif2 < 10. This is because M1/2 dominates 
the renormalization group behaviour. 

In this paper, we only address gravity mediated su- 
persymmetry breaking and do not consider other scenar- 
ios, such as gauge (GMSB) 78] or anomaly mediated 
(AMSB) supersymmetry breaking. Although, the 
low energy spectrum formulae we displayed in the pre- 
vious section are unchanged, the results for the bounds 
on the ^p-couplings or the LSP content change dramat- 
ically from one model to the other as we will see shortly. 
We hope that this paper serves as a basis to study the 
phenomenology of other SUSY breaking models. 

VII. RESULTS 

In the following numerical analysis, we use a version 
of SOFTSUSY [s^l which has been augmented with ^p- 
couplings. The beta functions for the .^p-MSSM cou- 
plings and masses contain the full one-loop and RFC 
contributions. The beta functions for the RFC MSSM 
couplings and masses also contain the two-loop pure 
RFC corrections. As discussed in Sect.V, small neutrino 
masses imply that the sneutrino vev's must be small. Al- 
though we derive their values from the minimization of 
the scalar potential, we neglect them in the calculation 
of sparticle masses. This is a good approximation, valid 
to 0{vi/MsuSY) ^ 1, when considering only the spec- 
trum of sparticles and not the small mixing induced by 
^p-couplings. We have checked that the error induced 
in the sparticle masses is much smaller than the current 
theoretical uncertainty in the RFC part of the calcula- 
tion mil Hi- The ^p contribution to the SM Yukawa 
couplings and fermion masses, however, is taken into ac- 
count as described in Sect. ]V\ Radiative electroweak 
symmetry breaking and the determination of sneutrino 
vev's follows the discussion in Sect. IIVI SOFTSUSY adds 
one-loop RFC threshold corrections to the sparticle and 
Higgs masses, and takes one-loop RFC threshold cor- 
rections into account when calculating the Yukawa and 
gauge couplings. For further details on the RFC part 
of the calculation consult Ref. |80j. Numerical results 
from the aumented version of the program SOFTSUSY, i.e. 
beta functions, neutrino masses, electroweak breaking, 
the mass spectrum, and bounds on the couplings etc 
have been carefully checked with an independent Fortran 
code 

We use the input parameters [s^ rrit — 174.3 GeV, 
aY^{Mz) = 0.1172 and mf,(m6)^ = 4.25 GeV, corre- 
sponding to mj^°'^ = 5.0 GeV at the 3-loop level. Other 
SM WS masses input are: m„(2 GeV) = 3.0 x 10"^ 
GeV, mc(mc) = 1.2 GeV, md{2 GeV) = 6.75 x 10~^ 
TOs(2 GeV) = 0.1175 GeV. The pole lepton masses are 
taken as nie = 5.11 x 10~* GeV, = 0.10566 GeV and 
rrir = 1.777 GeV. The Fermi constant Gp = 1.16637 x 
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10 ^ GeV ^, the fine structure constant a(0) = 
137.03599976 and Mz = 91.1876 GeV are used to de- 
termine the electroweak gauge couphngs. 

A. Bounds on Lepton-number Violating Couplings 

1. Procedure 

We first use the numerical analysis of the RGEs to 
set bounds upon the lepton-number violating couplings 
(Xijk, X'ijk) from the cosmological neutrino mass bound 
and requiring the absence of negative mass-squared 
scalars other than the Higgs and sneutrinos. (This does 
not refer to the physical mass and thus does not consti- 
tute a tachyon.) Neutrinos contribute to the hot dark 
matter and as such can free-stream out of smaller scale 
fluctuations during matter domination in the early uni- 
verse. This changes the shape of the matter power spec- 
trum and suppresses the amplitude of fluctuations. Com- 
bining the 2dFGRS data ^ together with the WMAP 
measurement 27] one can thus set a bound on the neu- 
trino mass at 95% C.L. 

^m^. <0.71eV. (99) 

i 

Scalar mass squared values can be driven negative 
during the RG evolution between the GUT- and the 
weak-scale, as happens to the Higgs in radiative elec- 
troweak symmetry breaking. But if any of the electrically 
charged or colour MSSM scalar fields develop negative 



mass squared values, QED or QCD would be broken, in 
conflict with observation. We therefore reject such values 
of A, A'. 

Neutrino mass and charge- and colour-breaking min- 
ima bounds depend not only upon the _^p-couplings, but 
also on the RPC SUSY breaking parameters. For a defi- 
nite quantitative analysis, we therefore take an example 
set of SUSY breaking parameters. We choose the SPSla 
mSUGRA point [s^ which has the following parameter 
values: Mo=100 GeV, M1/2 = 250 GeV, and trihnear 
couplings Ao = -100 GeV at Mx- tan/3(Afz) = 10 and 
/i > are also imposed. 

As stated in Sect. 13 a single non-zero ^p-coupling at 
Mx will generate through the coupled RGEs non-zero 
K,{Mz), D,{Mz) and (m^ ^XMz)- This is seen ex- 
plicitly in the RGEs in Eqs.' ||^[^IB31IbT6|) . where the 
anomalous dimension couples and well as the 
soft breaking sfermion masses, e.q. m^, with (m^, - ). 
Since the anomalous dimension 

7f; cx {YeAe + YdAd) , (100) 

Ki{Mz), Di{Mz), {m^^^j^ ){Mz) are also proportional to 

(YeAe+YdAd). Through k„ A, (m^^- J 9^ at the 
weak scale, we obtain non-zero sneutrino vev's, as can be 
seen from Eq. H53|) . This in turn gives us a non-zero neu- 
trino mass as seen in Eq. (|88|l . In order to estimate the 
resulting neutrino mass, we naively integrate the RGEs 
assuming constant parameters and insert our result into 
Eq. We obtain 



167rckGUT 



Vd 

167r2 



2 r 



In 



Mgut 
Mz 



El 



3A- ■ • {YD)jq + Xijq ■ {Ye)j 



11'^ 42 



B 

-Ml 



; tan/3) , (101) 



where / is a complicated dimensionless function of the 
SUSY parameters with typical values 0(10). A similar 
result was obtained some years ago by Nardi |2^. In 
Eq. IjlOlll . we explicitly see the dependence of the in- 
duced neutrino mass on the product of l^p- and Higgs- 
Yukawa couplings from Eq. I|100l) . Given a neutrino mass 
bound, e.g. Eq. (|99() . we can thus derive bounds on the 
^p-couplings. In the case where the down-like quark or 
the charged lepton mass matrix are diagonal, only the 
-^p-couplings A^j,j, or Ai^fc induce neutrino masses. Thus 
in the case of the LLi?-operators, since we do not include 
lepton mixing, we only obtain bounds on \ikk^ c-f- Ta- 
ble ^3 For the quarks we include the CKM-mixing and 
thus obtain bounds on all A', c.f. Table iHll 

Eq. H101|) works as an order of magnitude estimate. 
Setting aouT = 0.041, M1/2 = 250 GeV, tan/3 = 10, 
Yb = 0.18 and / = 10 and using the WMAP bound 



Eq. (jnni), we obtain 

^ ('3A^, • {Yn)jq + A„, • {YE)jq) < 2 • 10"^ . (102) 

With Yb ~ 0.18, we thus obtain the single bound 
A333 < 3 X 10"^. Full numerical integration shows that 
A333 < 6 X 10~®. Note that the only tan/3 dependence in 
Eq. H101|) is in the function /. 

Another interesting remark arises from Eq. H101|l : the 
higher the ultraviolet scale is (here denoted as Mqut) 
the larger the resulting neutrino mass and the stronger 
the bound on the A', A. Therefore, for the mSUGRA 
scenario, Mqut — 2 • 10"^^ GeV the bounds are stronger 
than for the GMSB model where Mqut must be taken 
at the intermediate energies 10^^ GeV. 

We also have to remark here on another independent 
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source for neutrino masses in the -^p-mSUGRA scenario 
coming from finite threshold effects involving squark- 
quark or slepton-lepton loops. The resulting neutrino 
masses are given by [gF 



2,2-K 



rr X! ^iki^jikml sin2(/)^ In 



'"I 



k.l 



d 

Kki^'jikmt sin 20^ In , 



k,l 



(103) 



with i^t) the lepton (down-quark) masses, (f)^{4''^) 
the slepton (squark) mixing angles and mj. (wj^ ) are 
the slepton (squark) mass eigenstates "88^. More details 
are found in Ref. (66> .87.] . Since the mixing in the first 
and second generation is negligible and also sleptons are 
almost degenerate the finite neutrino effects of Eq. 1] 103(1 
are not significant for the heaviest neutrino as compared 
to the ones induced from Eq. HlOlf) . For the third gener- 
ation we find 



^loop 



In- 



Ml/2 



3mi [ In 

TT \^ Mz 



(104) 



P 



The above estimate shows that bounds derived from 
Eg. HlOlfl are stronger than those derived from Eq. H103|) 
|89| . Thus the new bounds on the -^p-couplings presented 
in Table II are determined using the constraint Eq. (|99|l . 
the full solution to the one-loop RGEs and an accu- 
rate numerical diagonalisation of the ncutralino/neutrino 
mass matrix. 



coupled. Thus given one coupling at Mx in the weak ba- 
sis, we will in general generate an entire set at Mz (in the 
weak basis). In order to perform this computation, we 
must know the explicit form for the Higgs Yukawa ma- 
trices. However experimentally, all we know is the CKM 
matrix at the weak scale 



VcKM = UJ^Dl 



(108) 



as well as the diagonal matrices in the mass eigenstate 
basis. 

[mJ(jjj^g(Mz) = diag (rrid, TO^, TOf,)(Mz) , (109) 
[m«]diag(^^) ^ diag(m„, mc, mt)(Mz) . (110) 

For VcKM , we use the central values of the mixing angles 
in the "standard" parameterization detailed in Ref. j85( 

si2 = 0.2195, S23 = 0.039, S13 = 0.0031. (Ill) 

We neglect the CP-violating phase ^13 = 0. 

In order to perform the computation, we shall make 
the following simplifying assumptions. 

1. Due to the uncertainty concerning the neutrino 
masses and mixings we shall here assume that Y^; 
is diagonal in the weak current basis and thus 



(EL,R)ij — % 



(112) 



We shall return to the discussion of massive neutri- 
nos and their mixings in our framework in a future 
publication. 

2. We shall assume that Yo^jj are real and symmetric. 
Thus Ul = Ur and Dl = Dr. 



2. Quark Bases 

Before discussing our results, we must insert a dis- 
cussion on bases. In our initial parameter set at the 
GUT scale (c./. Eq. (jHII)), the -^p-coupHngs are given 
in the weak-current eigenstate basis. Similarly the Higgs 
Yukawa coupling matrices Yg, Y^i, Yjj, and the corre- 
sponding mass matrices are also given in this basis, i.e. 
in general they are not diagonal. The matrices are diag- 
onalized by rotating the left- and right-handed charged 
lepton and quark fields from the weak basis (") to the 
mass basis (™) 



(El.r)., (elj,) 
(Ul,r)„- 



3 1 



(105) 
(106) 
(107) 



In general, the rotation of the left-handed fields [e.g. 
Ul) is different from the right-handed fields (Ur). In 
the weak basis, due to the non-diagonal elements in 
Yb, Yu, Y[/, the RGEs for different ^p-couplings are 



3. When determining bounds below, we consider three 
extreme cases: (a) no-mixing, (b) the mixing is only 
in the down quark sector, (c) the mixing is only in 
the up-quark sector. This corresponds to 



(a) Dl,r = 1, 

{b) Dl,r = Vc_ft:Af, 

(c) Ul,r = VcKj\/, 



Ul,r = 


1, 


Ul,r = 


1, 


Dl,r = 


1 . 



(113) 



In these three scenarios, the mass matrices at the weak 
scale and in the weak current basis are then given by 



(a) md{Mz) 
m.u{Mz) 

(&) md(Afz) 
m„(Mz) 

(c) ma{Mz) 
niuiMz) 



["^dUiagiMz) , 



CKM 



■ [mrf]diag(^^^) • K 



T 

CKM J 



[m„]diag(A^2) . 
[^d]diag{Mz) , 
^CKM ■ [^u](\ia.g{Mz) 



^CKM ' 



(114) 
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Thus in each scenario, the matrices md{Mz), mu{Mz) 
are determined uniquely in terms of their eigenvalues and 
the CKM matrix. 

The Higgs Yukawa matrices y^D,u are proportional 
to the mass matrices. Therefore in each scenario of 
Eqs. H113I114(I the RGEs are fully determined. Given 
a set of -^p-couplings at Mx (of which we will here only 
choose one to be non-zero), we can then compute the ^p- 
couplings (including Ki) at the weak scale in the weak 
current basis. Given the full set of parameters at Mz we 
can diagonalize the neutrino/neutralino mass matrix in 
Eq. (|81|l and compute the neutrino mass. For a check this 
neutrino mass should be identical with the one derived 
in Eq. H86(l . We can then use the experimental bound 
on the neutrino mass, Eq. H99|l . to determine a bound on 
the -^p-coupling, in the weak current basis. 

For comparison with experiment we must rotate to 
the quark mass eigenstate bases in scenarios (6), (c) , 
Eq. (|113|l . To do this, we follow the procedure of Ref. |90l |. 
For scenario (6), with all the mixing in the down quark 
sector, we obtain the ^p-interactions for the superfields 
in the quark mass eigenbasis 



N^{VcKM)llDl - E,Uj 



+ 2'^ijki^CKM)mj{VcKM)nkUiDmDn- (115) 

Referring to Eq. (|115|l . we define the rotation of the cou- 
plings to the quark mass basis (denoted with a tilde) 



Kjmi'^CKM)mk, 



K'jk — ^imni'^CKM)mj{VcKM)nk- 



(116) 
(117) 



For scenario (c), with all mixing in the up-sector, and 
the superfields in the quark mass eigenstate basis, the 
superpotential terms are 



(&) 



D A' 



ijk 



N,D,-E,Ui{V, 



ckm)m 



+ -j^>^ijk{VcKM)liUlDjDk ■ 

This implies the rotation of _^p-couplings 

Kjk — Klki^CKM)jl 1 
Kjk = K']kiVcKM)il , 



(118) 

(119) 
(120) 



where in the first term we have taken the rotation of the 
EUD term. 

Another set of bounds applied on the -^p-couplings X'-j^ 
arises from the requirement of no sneutrino tachyons, i.e. 
we require the physical mass m? > 0. The resulting 
bound has been observed first by de Carlos and White 
and can be estimated as 



jk 



0.5M2,2-|- iM|cos2/3 



IStoq 



49M2/2 



12^2 



(121) 



For the SPSla benchmark scenario this bound sets all 
\[^j,{Mx) to be less than 0.13 in good agreement with 
the exact numerical solutions of the RGEs in Table II 
below. 



3. Discussion of the Bounds 

Table Hnl displavs the strongest upper bounds upon tri- 
linear A' couplings coming either from the neutrino mass 
constraint or the absence of tachyons at mSUGRA point 
SPSla as described in Sect. IVII A ll abovc. The different 
bounds coming from altering the quark mixing assump- 
tion are displayed. In each case, the upper bound at 
Mgut is shown in the weak eigenbasis, and the corre- 
sponding bound that is obtained when the couplings and 
masses of the MSSM are run down to Mz and rotated to 
the quark mass eigenbasis as in Eqs. H116I117I119I12()(I . 
Neglecting quark mixing we see that some of the bounds 
come from the absence of tachyons, and allow large 
couphngs of around 0.4 at Mz. However, for A-^^-, 
the diagonal components of produce a non-zero k 
through the RGEs, which in turn generates a neutrino 
mass. These bounds are much stronger and are of order 
O(10"3 - 10"^). It should be noted that the neutrino 
bounds are sensitive to the down quark mass inputs, be- 
cause the RGEs generate k proportional to Y^. When 
the CKM mixing is assumed to be in the up-quark sector, 
K23 ' Aii3 and A^gj acquire stronger bounds coming from 
neutrino masses because the larger up-quark Yukawa cou- 
plings in Y(7 also begin to mix the Y^ through the 
RGEs. When all down quarks are mixed at Mz, any X'^-^. 
produces n terms and therefore a non-zero neutrino mass. 
In this case, all of the bounds are strong: 0(10"'^ — 10^^). 

Table IIVI shows the equivalent bounds for the A pa- 
rameters. These bounds are not sensitive to assumptions 
about quark mixing because the RGE generation of k 
proceeds through the charged lepton Yukawa couplings, 
which we have assumed to be diagonal in the weak ba- 
sis at Mz. Changing this assumption should drastically 
change the presented results. We see that 3 of the 9 A 
couplings are not very strongly constrained; they are al- 
lowed to be 0(1). If the Y^; were strongly mixed, this 
would no longer be the case and the neutrino mass con- 
straint would provide stronger constraints, which we ex- 
pect to be at the level of 0(10"^) - 0(10"'^), similar to 
the 6 couplings that are constrained by neutrino masses 
in Table ITVI 

We may ask how much the bounds in Tables IIIIIIVI 
depend upon the supersymmetry breaking parameters. 
In order to investigate this issue, we scan over the pa- 
rameters of the no-scale mSUGRA a simple hy- 
persurface of mSUGRA parameter space where mo = 
A{) — 0. The remaining parameters (tan/3 and M1/2) 
are varied in Fig. ^ and the maximum possible value of 
log(A333(MG[/j-)) is displayed as the background colour, 
as referenced by the bar on the right hand side. The 
white region marked "no model" has tachyons for any 
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TABLE III: Upper bounds upon trilinear A' couplings for SPSla in the quark mass eigenbasis at the weak scale Mz and in 
the weak eigenbasis at the GUT scale Mgut- The quark mixing assumption is shown in the first row for each case. Input 
parameters are given in the text. A superscript of t, v denotes the fact that the strongest bound comes from the absence of 
tachyons or the neutrino mass constraint respectively. 



value of A333 and so is not valid. White contours of 
A(j33(TOaa;) = 10-^ lO'^ ^, lO""^-^ and IQ-^-^^ are shown 
from bottom to top respectively. The strongest bound 
comes from the neutrino mass constraint, and we see a 
variation of 2 orders of magnitude on the bound across 
the parameter space, the strongest bounds coming from 
the low Ml/2 region. The reader should note the M1/2 
dependence of the neutrino mass in the simple formula 
Eq. (|88|l . This strong variation of the neutrino bound is 
also apparent for the case of other A' couplings. Fig. [21 
shows the variation of the upper bound on A23i(Mg;7t) 
with no-scale mSUGRA parameter point. The strongest 
bound comes from the no tachyon constraint, and we see 
only a small variation of the bound across the param- 
eter space, the strongest bounds coming from the high 
tan/3 region, at low A/1/2 • (Recall the A/i/2 sensitiv- 
ity in Eq. (|121|) .l The behaviour of small variation in the 
tachyon bound with supersymmetry breaking parameters 



is replicated for other lepton-number violating couplings. 
The weak bound of k, 0.5 over much of the parameter 
space is dependent upon the no-charged lepton mixing 
at Mz assumption. 

It is instructive to compare the bounds derived here in 
a representative scenario of mSUGRA in Tables IIIII liy | 
with the 2ct bounds at Mz collected in Table 1 in Ref. 
for a rather generic R-parity violating scenario. For the 
comparison we choose the no mixing scenario, i.e. case 
(a) in Eqs. H113I114() and squark and slepton masses of 
order of 100 GeV in the latter. For the \[jf.LiQjDk cou- 
plings, we obtain here one order of magnitude improve- 
ment for A211, two orders of magnitude for A3ii,A'i22, 
three orders of magnitude for A'133, four orders of mag- 
nitude for A222JA322, five and up to six (!) orders of 
magnitude for A233 , A333 . The sneutrino tachyon con- 
straint of Eq. H121|) sets slightly stronger bounds on the 
couplings A323, A223,A232, Ai32,A33i. In the case of the 
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TABLE IV: Upper bounds upon trilinear A couplings for 
SPSla at the weak scale Mz and at the GUT scale Mqut- 
Input parameters are given in the text. A superscript of t, v 
denotes the fact that the strongest bound comes from the 
absence of tachyons or neutrino masses respectively. 




FIG. 1: Upper bound upon A333(Mg!7t) as a function of no- 
scale mSUGRA parameter point, assuming all quark mixing 
resides in the down sector at the weak scale. The background 
colour displays the bound as measured by the bar on the right 
hand side. Contours of iso-bound are also shown. In the top 
left-hand white region there is no tachyon-free model for any 
value of the coupling. 



-^p-couplings XijkLiLjEk we obtain two order of magni- 
tude stronger bounds than in Ref. j3| for the couplings: 
A122, A322, A133, A233. Sneutrino tachyons do not set bet- 
ter limits in this case. Comparison of the quark mixing 
cases (b) or (c) of Eqs. aii:^lll4|l derived in Table. Hill with 
the Table. IV of Ref. [l8j show similar orders of magni- 
tude, but stronger bounds for some of the couplings. 
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FIG. 2: Upper bound upon A23i(Mgi7t) as a function of no- 
scale mSUGRA parameter point. The background colour dis- 
plays the bound as measured by the bar on the right hand 
side. Contours of iso-bound are also shown. In the top left- 
hand white region there is no tachyon-free model for any value 
of the coupling. 



B. LSP Content in the No-Scale Model 

As outlined in the introduction, in ^p-mSUGRA the 
|?p-couplings can affect the weak-scale particle mass spec- 
trum via the RGEs. They can also affect the interpreta- 
tion of the resulting spectrum, since with the LSP is 
no longer stable, and thus no longer subject to cosmolog- 
ical constraints on stable relics. In the ^p-mSUGRA the 
LSP need not be electrically and colour neutral. Before 
discussing the ^p-case we briefly review the RPC case. 



1. The RPC Case 

To begin with, we perform the scan in the free pa- 
rameters Ml/2 and tan/3 in R-parity conserved no-scale 
mSUGRA. The LSP mass and contours of equal lightest- 
Higgs mass are displayed in Fig. |31 The background 
colour displays the LSP mass according to the scale on 
the right hand side of the plot. The region disallowed by 
tachyons is shown in black. In the bottom left-hand side 
of the plot is a white line which shows the boundary of the 
LSP identity. Below the line, the LSP is the lightest neu- 
tralino, whereas above it the LSP is a right-handed stau. 
A charged LSP is ruled out in the R-parity conserved 
scenario from cosmological constraints, and so the en- 
tire region above the white line is ruled out. This bound 
comes from limits on abundances of anomalously heavy 
isotopes ■ LEP2 |^ places a lower bound on the 
Standard Model Higgs mass of nih > 114.4 GeV. This can 
also be applied to the MSSM Higgs when sin(Q; — /3) « 1, 
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FIG. 3; No-scale supergravity in the R-parity conserved limit. 
Labeled constraints coming from tachyons are shown. The 
background colour displays the LSP mass, which can be read 
off from the bar on the right hand side. Dashed contours are 
contours of lightest Higgs mass. The white line delineates the 
labeled regions of f LSP and x? LSP. 



which is the case in all of our results. The theoretical 
uncertainty upon the lightest Higgs mass is estimated 
to be ±3 GeV [93, so we place a cautious lower bound 
on SDFTSUSY's prediction of 111 GeV. Even so, we see 
from Fig. 13 that there is no parameter space left with 
both a heavy enough Higgs and a neutral LSP. Thus no- 
scale supergravity is ruled out for the R-parity conserved 
MSSM. However, even a very tiny ^p-coupling will make 
the LSP unstable on cosmological time-scales and the 
neutral LSP constraint is then no longer applicable. For 
small couplings < 0(0.1), the spectrum can be approxi- 
mated by the R-parity conserved case, and so Fig.|31can 
still be used. We see that the entire region above the 
Higgs mass contour of 111 GeV would be allowed, for 
stau LSP masses above 96 GeV fgsj. 



2. The ftp-Case 

We now map out some parts of no-scale mSUGRA for 
-^1/2 — 500 GeV. Because we wish to show the effects of 
R-parity violation on the spectrum, we pick cases where 
the upper bound on the .^p-trilinear coupling is weak. 
This obviously occurs when the tachyon bound is the 
stronger of the two we have shown in Tables IHII and IIVI 
We display one A-type coupling (Fig^, one of type A' 
(FigEl) and one of type A" (FigEJ. 

Fig.Qlshows the variation of the nature of the LSP with 
tan/J and X23i{Mgut)- The case (b) of Eqs. H113I114|) is 
considered. For M1/2 = 500 GeV, as assumed here, we 
see from the equal Higgs mass contours, that the lower 
bound of 111 GeV on the lightest-Higgs mass does not 
pose a very severe constraint for tan/3 > 3. The LSP 




tan p 

FIG. 4: LSP content of no-scale mSUGRA for M1/2 = 500 
GeV, A231 non-zero at Mqut and weak-scale mixing entirely 
in the down quarks. The mass of the LSP is displayed in 
the background and corresponds to the bar on the right hand 
side. Regions ruled out by the presence of tachyons are in 
black. The white line delineates labeled regions of different 
LSP content (e for selectrons and r for staus). The dashed 
lines display contours of equal lightest Higgs mass. 



mass varies up to 190 GeV in the plane, but this value 
is a function of Mi/2. The diagonal white line separates 
regions of selectron LSP (above the white line) and stau 
LSP (below the white line). Note that there is an in- 
dependent (2a) bound for the coupling A231 from the 
known ratios Rr ~ r(r evD)/T[T — > ^lvv) correspond- 
ing to [Hl^: A23i(Mgut) < 0.046 x (mg^/lOO GeV). 
Comparing this bound with the nature of the LSP in 
Fig. 0] we observe that the scalar tau LSP is favoured for 
tan /? ~ 4 unless the above laboratory bound is evaded 
by taking M1/2 > 500 GeV. 

In Fig. we show the variation of the non-zero neu- 
trino mass in the tan/3 — A231 (Mg[/t) plane. Neutrino 
masses provide the upper bound upon A231 for mixing in 
the up-quark sector [case (c) in Eqs. (|113lll4|l ]. as as- 
sumed here. For larger values of X^'^ ~ 0.15, neutrino 
masses of 0(0.1 eV) are possible. In this case, above 
the white line, the LSP is a tau sneutrino, and below 
it the LSP is the stau. The laboratory bound for the 
coupling A23i(Mgut) reads jl^ : A23i(Mgut) < 
0.057 X (mg /lOO GeV) and since we find that for the in- 
puts of Fig.tSlthe bottom squark mass is about 1.2 TeV, 
the laboratory bound is evaded: the stronger bound on 
A23]^ comes from the sneutrino tachyon as is shown in the 
upper half of Fig. [S] 

Finally, we investigate the case of baryon number viola- 
tion. The case (b) of Eqs. H113I114|I is considered. Fig. El 
shows how the no-scale mSUGRA LSP mass varies with 
tan/3 and X';^23(Mgut)- There is little variation with the 
^p-coupling, contrary to the lightest Higgs mass, which 
is displayed in the form of contours. The previous bound 
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FIG. 5: LSP content of no-scale mSUGRA for Afi/a = 500 
GeV, A231 non-zero at Mqut and weak-scale mixing entirely 
in the up quarks. The logarithm of the mass of the heaviest 
neutrino is displayed in the background and corresponds to 
the bar on the right-hand side. Regions ruled out by the 
presence of tachyons are in black. The white line delineates 
labeled regions of different LSP content. The dashed lines 
display contours of equal lightest Higgs mass. 




FIG. 6: LSP content of no-scale mSUGRA for M1/2 = 500 
GeV, A323 non-zero at Mqut and weak-scale quark mixing 
in the down sector. The mass of the LSP is displayed in the 
background and corresponds to the bar on the right hand 
side. Regions ruled out by the presence of tachyons are in 
black. There is a stau LSP throughout all of the parameter 
space. The dashed lines display contours of equal lightest 
Higgs mass. 



on A323 (see Table IV of bound (ij) apart from the the- 
oretical perturbativity bound comes from the leptonic Z- 
width ratio and is A323(Mgut) < 0.015 for quark mixing 
solely in the down-quark sector, and with a little varia- 
tion from A/1/2. We observe from Fig. |H1 that the stau is 
again the LSP. 

We have exhibited, in Figs. viable regions of 

MSSM parameter space where the LSP is the selectron, 
the stau or the stau sneutrino. Different LSP content 
drastically alters the collider signatures of the models. 
The analysis above showed a preference to the stau be- 
ing the LSP. We discuss this in some more detail in 
Sect. rvml below. 



C. Sneutrino- Antisneutrino Mixing with Stau LSP 

Models which violate lepton number by two units 
(AL = 2) and generate neutrino masses, also result in 
a mass splitting of scalar neutrinos and anti-neutrinos of 
the same flavour usually referred in the literature as sneu- 
trino anti-sneutrino mixing |95l lo^ lOTj . If the sneutrino 
mass difference Amp = moj^ — mj}_ , is large and the 
sneutrino branching ratio into a charged lepton is exper- 
imentally significant, then a like sign-dilepton signal in 
6+ e~ — *■ i>- i>+ with i> ^ l~ + X could be observed [95l| . 
Like the B-meson mass splitting , the observability of the 
sneutrino mixing effects depend on the ratio 



Amp 



(122) 



where Fp is the total sneutrino decay rate. As we have 
already seen from Figs. 14161 in the no-scale scenario the 
stau, f , is the LSP when the -^p-couplings are small. In 
this (approximately RFC) case the specific flavour i — 
(e, fi) sneutrino decays, via charginos and neutralinos 
into i^e £~ f"*" i^r and C'g I'g r^. In this case, the 
probability of tagging a like-sign dilepton in the process 

1)1 i-T+i^risv{e^e^) = v {£+£+) +vie~e-) with ^ 



2(1 



(123) 



We investigate below the magnitude of this probability in 
the no-scale model with Mi/ 2 — 500 GeV and tan /3 = 20 
and with one dominant _^p-coupling A'i22(AfGUT) = 7.5 x 
10~^. Furthermore we consider no-quark mixing in de- 
termining the relevant bounds from the neutrino masses. 
In this model the stau is the LSP. We first calculate the 
sneutrino mass squared difference 



Amp — 



2m r, 



2mp 



(124) 



where mp is the average mass of mp^ . The sneutrino 
mass difference has been calculated in Rcf. 98] in a 
general basis independent manner. With our choice 
\'i22 generate at the electroweak scale the non-zero 
|?p-parameter set: wi, ki, Di, (m^ - ). The other l^p- 

parameters remain zero |99|. This simplifies our calcu- 
lation for the sneutrino mass splitting, since we can use 
the case of one sneutrino generation (the other two decou- 
ple from the mass matrices Eqs. (|59l61f) ). The sneutrino 
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mass splitting reads |98| : 

2 M| sin^ P sin^ ,5 



{Mf^„-ml){Aq„^ml){Ml„~ml) ' 



with 



cos (5 — ±- 



(125) 



(126) 



Notice that Eq. (|125|l does not depend on the super- 
potential parameters in contrast to the neutrino mass 
in E g. It is helpful to see the numerical values 

|lO0j| for the parameters at the electroweak scale start- 
ing from the no-scale model defined by: M1/2 — 500 
GeV, tan/3 = 20 and A'i22(Mgut) = 7.5 x 10"^ We ob- 
tain: B{Mz) = 33 238 GeV^ = 357 GeV, M^o = 91 
GeV, Mho = 816 GeV, M^o = 816^GeV, Vd{Mz) = 8.7 
GeV, vi{Mz) = -0.0012 GeV, Di ^ -0.74 GcV^ 
and (m^^-^) = 2.5 GeV^. Applying these values to 

Eqs. H125I126|I we obtain sin^ (5 = 1.3 x 10"^ and Amp = 
2.5 eV. The sneutrino mass splitting is of the same order 
as the neutrino mass obtained from Eq. (|86|l . since for 
fi{Mz) = 817 GeV and ki{Mz) = 3.5 x 10"^ GeV we 
have = 1.2 eV 101] . 

In order to calculate the probability P(£^£^) we still 
need the total sneutrino decay rate and the branching 
ratio B{h'e £~f'^i'r)- In the above scenario the right 
handed selectron of the third generation (we call it stau 
here although it is in fact an admixture of the three 
charged sleptons with the charged Higgs boson states) 
is the LSP with a mass rrif = 162 GeV. The rates for 
the chargino and neutralino mediated sneutrino decays 
(which we assume to be the dominant ones) are |95j| : 



with 



g2mlm'^ tan^ f3f^+ {rrif/ml) 
15367r3(M2, sin 2/3 - Mz/i)^ ' 



g^m|/^o(m?/m|) 
30727r3M4 



(127) 



f^+ix) = (1-.t) (l + 10a; + a;^) + 6x(l + a;)lna;, 



Vlx^ Ini 



(128) 



In the no-scale model under consideration we obtain: 
Mi{Mz) = 206 GeV, AhiMz) = 411 GeV with the 
gauge couplings giMz) = 0.3574 and g2{Mz) = 0.6525. 
Thus from Eq. (fn7|) we obtain: T[vt i-f+Vr) = 210 
eVandr(i>^ v^t^t^) = 1.2xl05eV. Sox,. = 2x10"^ 
and Bivi l'^ f + Vr) — 1.7 x 10^"^. We conclude that in 
this numerical example the probability for like sign dilep- 
tons, Eq. fl^ . is: Vii^l'^) = 6 x 10-l^ far too smaU 
to be observable. Of course this result depends on the 
parameter space and the probability 'P(^^f^) is bigger 
for smaller values of M-yj^ and larger tan/3 values (see 



Eq. H127|) . However, if we take into account the current 

experimental data, Afi/2 ~ 200, then V{l^t'^) ~ 10"^ 
We obtain similar results for the other _^p-couplings. 

The above benchmark computation can be helpful to 
the reader in order understand the typical magnitude of 
the parameters we are dealing with in this paper. 



VIII. STAU-LSP PHENOMENOLOGY 

As discussed in Sect.U in the case of ^p, the LSP need 
not be the lightest neutralino, Xi- In the previous sec- 
tion we have investigated the nature of the LSP in the 
mSUGRA scenario and have found regions in parameter 
space with different LSP's. In Fig. ^ we have a selec- 
tron or stau LSP, in Fig. [S] we have found a tau sneutrino 
or a stau and in Fig. we have found a stau LSP. The 
bounds in Table IIIII imply that if there is any apprecia- 
ble CKM mixing in the down-quark sector at the weak 
scale, must be very small. We also see some strict 
bounds upon the A^fe in Table llVl If the ^p-couplings are 
very small, the spectrum has negligible perturbation from 
the R-parity conserved case, the LSP content of which is 
displayed in Fig. |3| The allowed parameter space with 
TO/iO > 111 GeV in Fig. El then leads to a stau LSP. Thus 
we see a preference for a stau LSP in many no-scale R- 
parity violating scenarios. 

In the RPC-MSSM, the collider phenomenology relies 
crucially on the Xi-LSP, with all produced sparticles de- 
caying in the detector to x5 plus other i?p-even particles. 
This results in missing transverse energy as a typical sig- 
nature for all production processes. In the -^p-MSSM the 
RGEs and thus the spectrum, is altered. This changes 
the decay chains. Since typically all decay chains end in 
the LSP, the nature of the LSP is essential in determining 
the supersymmetric signatures. A detailed investigation 
is beyond the scope of this paper. We shall here focus on 
a classification of the signatures for the main production 
processes in the case of a stau LSP. 



1. Stau Decays 

The following discussion of the stau-LSP is somewhat 
analogous to the discussion in Ref. j34j for the x5"LSP. 
In determining the final state signature it is important 
to know how the stau-LSP decays. We shall assume that 
there is a hierarchy among the ^p-coupling constants 
with one dominant coupling, similar to the SM Yukawa 
couplings in the mass eigenstate basis. We furthermore 
assume the mixing due to Ki is small as seen in the previ- 
ous sections of this paper. Then there are two important 
distinct cases. 

1. The stau couples to the dominant opera- 
tor. The dominant operator is in the set 
{Le,^iLrEe^^^T, L^Lf^Er, LrQiDj}. In this case, 
the stau simply decays via the two-body mode. 
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For the dominant operator LrQiDi for example 
we then obtain 



r(f- 



d) 



167r 



(129) 



where Nc — 3 is the number of cofours. The 
com plete list of ^p, two-body decays is given in 
Ref. |l02j | . For a recent treatment of two-body stau 
decays also see [s^l ■ For the above two-body decay 
mode the decay length is given by 



CTf = 3.3 10""m 



10- 



A 



311 



100 GeV \ 
Mf J 



(130) 



which in an experiment must be multiplied by the 
relevant Lorentz boost factor 7^ of the stau. Only 
for very small coupling (A' ^ 10~^) is the decay 
length relevant. 

2. The stau doesn't couple to the dominant 
operator. The dominant operator is in the set 



{LeLf_tEe,f_t, Le^f_,QtDj, UiDjDk}. In this case the 
f decays via a four-body mode. For the operator 
LfiQiDi there are four decay modes via the 
neutralino 



,+ 



Vn. + d + d 



u + d 
u + d 
d + d 



and three decay modes via the chargino 



'^T + (Xi )* -> l^T 




(131) 



(132) 



As an example we here compute the decay f — > 
T~ fi~ud. The details of the computation, in par- 
ticular the four-body phase space are given in Ap- 
pendix^ The result is 



J 



T{r- ^ r-^,+ud) = MUK\' 



— — M~ 

237r5M2m4 ^ ' 



(133) 



where K 



1/(720 X 2^) = 1/23040. ar,d, 



are neutralino coupling constants given in the ap- 
pendix, is the neutralino mass and m is the 
universal scalar fermion mass. We have assumed 
massless final state particles and neglected the mo- 



menta compared to Af^, rh. In the last step we have 
set the couplings Or^d — b^,u = g, the weak coupling 
constant. If the four-body decay is the dominant 
decay mode, the decay length can be estimated as 



J 



CTf = 6.2 IQ- 



10 



-3 



V 100 GeV J V 100 GeV J \ Mf 



^ flOO GcYV 



(134) 



r 



For reasonable supersymmetric masses and cou- 
plings this could lead to detached vertices in the 
detector. This is a very promising signature for the 
stau-LSP. 

If the two-body decay is allowed, i.e. the relevant cou- 
pling is not suppressed, it usually dominate over the four- 
body decay. In order to estimate the required hierarchy 
of couplings for the four-body decay to be relevant we 
consider the ratio 



r4(7 



T2{f — > ud) 



O 



3ij 



> 1. (135) 



Assuming the sparticle masses are roughly equal, this 

corresponds to A'2ii/A'3.y ^ 0(10^) for the 4-body decay 
mode to dominate over the 2-body one. If, for example, 
= rh = 2M?, we obtain A'2ii/A'3y ~ 0(10^), which 
is not an unreasonable hierarchy between generations. 



2. Collider Signatures 

At a collider, the main supersymmetric pair production 
processes are 



gg, qq, e+r, X°X°: 



Xi Xj 



Xi Xj • 



(136) 
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Here we investigate the possible signatures for these pro- 
cesses in the case of a stau LSP. In order to determine the 
final state within the detector, we must know the decay 
patterns of the particles. This strongly depends on the 
supersymmetric spectrum and thus upon which point in 
SUSY breaking parameter space is being studied. For 
this first study, we shall assume the mass ordering 



and the final state signature will be 



Trig > rriq > rriji > Tn^± > ni^o > rrif 



(137) 



which we typically obtain (with or without ^p) within 
mSUGRA. If there are no near-degenerate particles, a 
produced supersymmetric particle will dominantly cas- 
cade in two-particle decays down the mass chain (jl37ll . 
We display this decay chain in Fig.[7| We have added at 
the end both two- (in red) and four-particle (in blue) stau 
decays. Final state quarks are denoted by "j" to indicate 
a jet. We can use this decay chain to determine a quali- 



I X 





1,2 




FIG. 7: Possible dominant links in a sparticle decay chain 
with a stau LSP and R-parity violation. Two-body decay 
modes of the LSP are shown in red, and 4-body in blue. 



tative picture of the possible final state signatures. Note 
that due to the strict bounds on the -^p-couplings which 
we have determined we only expect these to be relevant 
in the stau-LSP decay. Furthermore, in determining sig- 
natures we shall assume that either the two-body or the 
four-body stau decay dominates. We do not consider the 
case of comparable partial decay widths. 

At the Tevatron and LHC the largest production cross 
sections are for gluinos and squarks. If we consider for ex- 
ample QBtlR production, then the dominant decay mode 
for the squark is 



QRqR 



6j 
6j 
6j 
6j 



an 



-2(T^ 
2(r+ 



for 
for 
for 
for 



■ JJ 

• Ti 
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Here, any charge combination for the leptons £ is allowed, 
since they result from the decay of a virtual (Majorana) 
neutralino (c./. App. D). This can give us a like-sign di- 
Icpton signature. Otherwise, we see that we have a large 
number of jets in the final state independent of the decay 
mode of the stau. (This would be reduced for dominant 
operators LLE.) This makes it more difficult to observe 
isolated high px charged leptons. We can obtain missing 
transverse momentum from the final state neutrinos but 
it will be extremely diluted due to the many body-decays. 
The most promising signature are like- sign dileptons to- 
gether with the direct detection of t's |l03l| . which is of 
course difficult. 

For qigi-production, we expect a larger liklihood for 
the cascade decay through the heavier neutralinos and 
also throu gh t he charginos. This can lead to a tri-lepton 
signature |l04j which can be extended by the additional 
t's. This requires a detailed analysis but we expect this 
to be more promising than the ququ outlined above. 

The gluino decays via the squarks adding an extra jet 
to the final state. In this case it might be more promis- 
ing to consider non-dominant decay modes, including a 
possible direct -^p-decay of the neutralino. An estimate 
of the relative rates for a pure wino neutralino is 



r(x;^M + 2j) 3A'^ 



3-10" 



(140) 



for A' < 10~^, and where we have neglected the stau 
mass. This is hopeless, unless the neutralino and the 
stau are nearly degenerate. 

At the Tevatron and LHC the pair production of slep- 
tons is about two to three orders of magnitude lower 
than the production of squarks or gluinos, for equal mass. 
However, we expect the mass to be lower (c./. Ea. (|137|l L 
and also the signal cleaner. At a future linear collider 
e^e~ facility this is typically an ideal mode for searches 
or the measurement of MSSM parameters. As we can see 
from the decay chain in Fig. d the slepton dominantly 
decays as 



We then obtain the final-state signatures 



(141) 



(138) 



4j + « + T-T+ 

4j + 2(«) + 2(T+r-) 
4j + M + W + 2{t+t 
4j + « + i/£ + 2(t+t- 



for f 

for f 

for f 

for f 



■ JJ 



(142) 
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In the second case the sign of the charge of the two lep- 
tons from the stau decays is arbitrary due to the interme- 
diate (Majorana) neutrahno. Thus we can have hke-sign 
tri-leptons, which is a very promising signature. 

Similarly, using the results from Fig. [TJ we expect the 
as dominant signatures for ncutralino pair production 




(143) 



depending on the decay of the stau-LSP decay which in 
turn depends on the dominant ^p-coupling. For chargino 
pair production we have 



Xi Xi 



4j + l/rVr 

4j + U+l/rl^T + T+T- 
4j + + VrVT + r+T" 



for f ^ jj 
for f —t r^jj 
for f TV}\ ' 
for f t[v,1)\\ 

(144) 

assuming the chargino decays directly to the stau LSP. If 
we produce the heavier electroweak gauginos we can cas- 
cade decay through the lighter gauginos producing more 
charged leptons. For the neutralino we have promising 
multi-lepton sigantures, whereas for the chargino we ex- 
pect a significant amount of missing ■ 

In summary, as promising signatures in the case of the 
stau LSP we have 

1. A detached vertex from the long lived stau, partic- 
ularly in the case of the four-body stau decay. 

2. Multilepton final states. 

3. Multi-tau final states, requiring efficient tau tag- 
ging. 

The four-body decay of the stau results in more final state 
leptons than the two-body decay and is thus possibly 
more promising. 



IX. SUMMARY AND CONCLUSIONS 

We have investigated for the first time the general 
MSSM in the context of mSUGRA. We have studied in 
some detail the origin of lepton-number violation and 
have found that with respect to the dimension-five oper- 
ators baryon-parity is preferred over R-parity. We have 
then shown that in a wide class of models both and Di 
are zero after supersymmetry breaking at the unification 
scale. We have taken this as our boundary conditions at 
Mx in order to investigate the resulting model in consid- 
erable detail. 

In order to embed the model within the unification 
picture we have computed the full set of renormalza- 
tion group equations in the appendices. We have used 



two methods, including a novel method of Jones et al, 
which is particularly conducive to the numerical imple- 
mentation. We then developed an iterative algorithm 
which solves the RGEs, minimises the potentail of the 
five neutral, scalar, CP-even fields, while implementing 
weak-scale Yukawa and gauge boundary conditions. The 
algorithm is stable and has been checked by an indepen- 
dent program. This is one of the main technical advances 
in this paper. Given the minimum, we determined the 
complete supersymmetric spectrum, including also the 
mass of the heaviest neutrino. 

We have then shown that the ^p-couplings in this 
model are severely constrained by the upper bound on 
the neutrino masses, as summarized in Tables IIIII and 
El Thus when embedding the ^p-MSSM in mSUGRA 
the neutrino mass bound is the strictest and most univer- 
sal, i.e. applies to all lepton number violating couplings. 
This is one of the main results of this paper. 

We have then looked in detail at the nature of the 
LSP. We have found solutions with a selectron, tau sneu- 
trino and stau LSP besides the usual neutralino LSP, 
with the stau most favoured in the noscale mSUGRA 
model. This significantly affects collider phenomenol- 
ogy. We present a first discussion of this broad topic 
in Sect. IVIIII We have also studied the phenomenology 
of sneutrino-antisneutrino mixing in this model, but do 
not expect any significant effect. 

We conclude that the -^p-MSSM is as viable as the 
RPC-MSSM. As we show, it considerably differs both 
conceptually and phenomenologically from the RPC. The 
intimate connection with neutrino masses is an outstand- 
ing feature which we shall discuss in more detail in a 
forthcoming publication. 
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APPENDIX A: NOTATION AND ANOMALOUS DIMENSIONS 



The chiral superfields of the i?p-MSSM and the ^p-MSSM have the following Gsm = SU{3)c x SU{2)l x U{1)y 
quantum numbers 



L: (1,2,-i), E: (1,1,1), Q:(3,2,i), [7: (3,1,^), 
D : (3,1,-i), H,: (1,2,-^), : (1,2,^). 

The -^p-MSSM superpotential is then given by 

W = tab [(YshL^H'lEj + [YDhQTH\Dr. + [Yu)^JQT H^2U3A ~ fab [A*^i-ff2 + «;*i"^2] 



(Al) 



-{AEk)i.jL°;L''jEk + {Ajyk)ijL'}Q^''Dkx 



(A2) 



We denote an SU{3) colour index of the fundamental representation by x,y,z = 1,2,3. The SU{2)l fundamental 
representation indices are denoted by a, 5, c = 1, 2 and the generation indices by i,j,k — 1, 2, 3. We have introduced 
the twelve 3x3 matrices 



Y_E, Yd, Yij, A^k, Ajjk, Ajji, 
for all the Yukawa couplings. This implies the following conventions in the Martin and Vaughn notation 



Y 



Lf L ■ El- 



Y 



Y 



— {■^E>')ij^ab — ^{■^E>')ji^ab — Kjk ^ab, 



i — 1 —1 — txyz [-li-U^ )jk ^ ^xyz y-i^U' )kj — ^xyz-^ijk: 

The soft SUSY breaking Lagrangian is given by, 

BHiH2 + D,L,H2 +H.C 



E{m^^)E^ + D{mf^^)D + U{m^^)U 



+ {h.Ek)ijLiLjEk + {h.Dk)ijLiQjDk + {\iu^)]kUiDjDk+ H.c 

where we have introduced the soft SUSY breaking trilinear couplings 

h_B, h/), hjy, h.Ek, h^ifc, h[/i. 



defined analogously as the Yukawa couplings in Ij A4IA6P . 

In general the one-loop renormalization group equations for the Yukawa couplings are given by [T^ 



dt 



167r 



7*= 

2 IP 



+ {k ^ i) + {k ^ j) . 



and the anomalous dimensions are 



\Y,„Y'P''-25lY,alCa{i).. 



(A3) 



(A4) 
(A5) 
(A6) 



(A7) 



(AS) 



(A9) 



(AlO) 



We have denoted by Ca( f) the quadratic Casimir of the representation / of the gauge group Ga- For details see the 
Appendix A of Ref. |l4| . All equations in this section are valid in the DR renormalisation scheme. 
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The one-loop anomalous dimensions are given by [ij, Il05| : 

7^; = [YEYl)^^ + {XE.A^E^h+-i{^D.A^D^h-5]{^9l + lgl). (AH) 

7|; = 2{Y%YE)^+TT{\E^K^^^)~5]{^-gl), (A12) 

7^; = (Yz^yI,)^^. + (Yc/Y[,)^^. + (Al,,Az,0.. - -^^l^^s? + ^52^ + ^5^) , (A13) 

7^?; = 2 (y^Yz,)^.^ + 2Tr(A^,A^0 + 2(A^,At,,),, - 5]{^gl + ^53^) , (A14) 

7^;; = 2 (y^,Y^)^.^ + Tr(A^.At,,) - 5]{^gl + ^-gl) , (A15) 

7^J = Tr (sYzjYl, + Y^Y^,) - + , (A16) 

7t = 3Tr(Y^Yt,)-(Ag? + |g|), (A17) 

^ 1h\ =-^i^D.^D)^,~{^*E.^E)^,■ (Al8) 



Note that here, Hi 2,L,Q represent the fields i?"2' where a is the index of the fundamental representation 

of SU{2) {i.e. no factors of eat are factored). The /3-functions for the Yukawa couplings 1^ and for the bilinear 
superpotential couplings are combinations of the above anomalous dimensions I)A11IA18(I . The two loop anomalous 
dimensions in the ^p-MSSM can be found in We present the one-loop beta functions for the superpotential 



couplings and masses for completeness. 

The RGEs for the Yukawa couplings (including full family dependence) are given by 

16n^^^{YEh = (yEhklEl + (JEhlSl - (As. )k^l^: + (YE)knt ' (A19) 

WTT^liYoh, = (YDhklDl + (ynhnHl - i^D. )k^7f: + (YD)kjlg , (A20) 

16n^±(Yuh = i'yu)rklu: + (YuhiSl + (yu)k,-f2: , (A21) 

i67r2A(A^.)..- - i^E'hii: + i^E^hnt: + {YeUih[ - i^E'^hit: - (YeUi^^ , (A22) 

167r2^(A^0.j- = (A^0.j-7Dr + (ADO.aS+(Az,.),,7L;-(Yi5),fe7i^;, (A23) 

16^2^(Ac,0,fe = {Au^hi^;; + {Au.)iklD: + i^u')jkl^: ■ (A24) 
Here t = ln(Q), and Q is the renormalization scale. The RGEs for the bilinear terms are 

16^2|/. = f^{lSl+lSl} + ^'lL:, (A25) 

16^2 rf. - ^'1^: + n^lll + ^'^H\ ■ (A26) 



APPENDIX B: A METHOD TO DERIVE THE SOFT SUSY BREAKING RGES 

A straightforward way to derive the RGEs for the soft SUSY breaking couplings and masses is by a direct use of 
the explicit formulae at 1-loop given in ^16j. This is a somewhat tedious job. A very elegant method which is also 
very helpful for numerical calculations is the one described in Ref . _15J ■ AH the soft SUS Y R GEs can be derived from 
the anomalous dimensions IjAlllAlSp by the action of an operator which is given below |106| . The method works not 
only at one loop but it has been proven to all orders in perturbation theory 1 151 . In principle one could apply the 
operators (B1-B5) below to the 2-loop anomalous dimensions derived in Ref. and write down the full two loop 
coupled RGEs in the most general case. However, here we restrict ourselves to the one-loop case. In particular the 
soft /3-functions for the bilinear b*^ , trilinear h*-''^' and scalar masses (m^)* soft SUSY breaking terms can be read 
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from 



IGtt 



dt 



dt 

d(m^) 



7;'b^'+7/b''-2(7i)l/i^'-2(7i)fM 
7,*h^'=' + 7/h*'^'' + 7f h^'*' - 2(7i)jYJ''=' - 2(7i)^Y*'=' - 2(7i)fY^''' , 



2\i 
j 



dt 



200* + 2MM*gi— + Y,,„„— 

dgi oYimn 



rlrnn 



where 



(71 



o 



(Bl) 
(B2) 
(B3) 

(B4) 
(B5) 



and repeated indices are summed over. At one loop the last term Xa ,in Eq. (|B3p is not relevant. Its (scheme 
dependent form) is given for example in the last reference of Ref. [I3 (see their Eq. (2.11)). The RGEs (|B1IB3|) are 
valid as long as we do not eliminate the U(l) Fayet-Iliopoulos (FI) D-term. The RGE running of the Fl-term can 
then be written independently. It is known that for universal boundary conditions this term is not renormalized down 
to low energies and we do not discuss its RGE here. On the other hand if we eliminate the FI D-term by using its 
equation of motion then this renormalization gives rise to additional contributions proportional to the U(l) gauge 
coupling (see the 5-tcrm in the RGEs for the soft SUSY breaking masses in the Appendix Q. Now from Eq. IjBip 
the RGEs for the bilinear soft SUSY breaking masses in the ^p-MSSM are 



167r- 



167r^ 



,_dB 

m 

dt 



B 



B1h\ 



{ll)H\ + ill)ll 



(7i)l:«' + (7i)h^« 



H2 .J 



2^(71)^; 



(B6) 
(B7) 



The RGEs for the trilinear soft SUSY breaking masses in the _^j,-MSSM can be read from Eq. HB2p 



dt 



16n- 



,dih 



Djik 



dt 



7^; {\^E)ik + 1h[ (^eU + ill (hi^-" )^l + Te: (hij)« 
2(7i)L;(YB)ifc - 2(7i)fj(Yi;),fc - 2(7i)f/(AsO«i " 2(7i)|'= (Y^),, 



7^; (hc)ife + il\ (hc),fc - 7^^ (hz,. );. + IdI k^D)a 



2(7i)^:(Yz5)/fc - 2[n^)l\{Yu\u + 2(7i)f;(A^.);, - 2(7i)^f (Y^),, 



IGtt 



2d(h[/)ifc 
~dt 



7^;(hc,),fe+7^^^(hy)a.+7[;r(h[/).i 
2(7i)g:(Yt/),fc ~ 2(7i)g^(Yy),fc - 2(71)^;^ (Yy),, 



(B8) 



(B9) 



(BIO) 



167r 



,d(hBfc) 



n _ Li 



dt 



IlI (hB^- )l3 - 1h\ (h£;)jfc + Tl, {^E'')iI + 1h\ {^E)^k + 7^f (hs' )^j 

2(71)^; (A^.),, + 2(7i)^^(Ys),fc - 2(7i)^;(A^.),i - 2(71)^^ (Y^),,- - 2(71)^ (A^,)..- ,(B11) 



167r 



2d( hjfc)a 

dt 



167r 



2 (i(h^.)-,-fc 
dt 



2(7i)^;(A^.)0- + 2(7i)^,(Yz5),fc - 2(7i)g(A^.),, - 2(71)^^ (A^Ou 
7^;; (h^;, ),, + 7^/ {\yu.)ik + (hu, ),i 



2(7i)^;(AyO.-fc - 2(7i)^^(A^0;fc - 2(7i)^f (Ac,0,7 • 



(B12) 



(B13) 
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The RGBs for the soft SUSY breaking masses in the ^p-MSSM can be obtained from Eq. HB3|I 
^Qy±^^^^^,d(rn^ = 4(ht,h^),. + 2Tr(h^.ht,,)-5,,(^5?|Mip 

+ 2{Y\,Y^),, + Tr(Al.A^,) + 2(YtjY£),, + Tr(AB.A;,) , (B14) 



+ (YeYI\j + (Ab,A^,)„- + 3(Az),A|,,)y , (B15) 



- 3(Al3,Yc),, - (A^,Yb),, - 3(A^,Yc),, - {K*^,YeU ^ (B16) 



(Yt;Yt,),, + (Al,, 
+ 2(Y|,Y^),, + 2 Tr(A^, A^O + 2(A^,aJ;,),, 

+ 2(Y|,Yi3)j, + 2 Tr(A^,ApO + 2(At;, A^,, , (B18) 



,^(n^^^^^,<^ . 4(ht,h.)..+2Tr(V.ht,J-.Jg,?|M,P + f, 



+ 2(Y^Yt;),, + Tr(A^.At,^) + 2(Yt,Yt;),, + Tr(At;.A^^,) , (B19) 

16^'^^ = Tr(6hchl, + 2hj5h^) - f -g?|Mi|2 + (Sgl\M^f 

+ 3 Tr(Y^Y|,) + Tr(YEYtj) + 3 Tr(Y,,Yt,) + Tr(Y£;Yt,) , (B20) 

+ 3 Tr(Y^Yt,) + 3 Tr(Yc;Yt,) , (B21) 

where from Eq. ljB4|) we have 

(7i)l;. = - (hijYl,)^^. - (hi,,A^,)., - 3(hz,,Al,,)y - b]{^^M^gl + \M^gl) , (B22) 

(7i)b; = -2(Y;hB)^^-Tr(hs.At,^.)-<5;.(^Afig?), (B23) 

(7i)g: = - (IidYU _ . - (h^Yt,) _ ^ - (Al,,h,,,)r» - ^){^r.^^9\ + ^Af^s' + ^^353) , (B24) 



(7i):C; = -2 (yLHd) .. - 2Tr(At,,.hoO - 2(hc;, A^,,),, - b){^M^gl + ^Mg^I) , (B25) 
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Tr (shoY^^ + h^Y^) - i^M.gj + hhgl) , 



-3Tr [huYl) - i^M.gl + hl^gl) 



(B26) 

(B27) 

(B28) 
(B29) 



and from Eq. ||B 



{YE)ik 


= (yE)ik{mi%i ^ 


f (Y£;),fcm^^ + (ABk),i(m^^^^) + 




)ii{m.f~'^)ik , 


(B30) 


{YD)ik 


= (Y£,)ifc(mQ^),i 


- (Apfc)H(m^^^^) + iYD)^km'H, + 


(Yi5),;(mo^)ifc , 


(B31) 


(yu)ik 




+ {Yu)tkm\^ + (Yc/)j((mu^);fc , 






(B32) 




= (A£;fc);j(mL^).a 


-(Ys),fe(m|^^^) + (As.).KmL') 


Jl + 


(YEhkiml^H,) + i^E')z3{rn^%k , 


(B33) 




= (A^fc)(j(mL^).a 


- (Yi3)j7c(m|^^^) + (Ack),,(mQ^ 




(A^O'u(mD^)ifc , 


(B34) 




= (Ac/Oife(mu^)i» 


+ (A[/Oifc(mf,^)ij + (Ac/Oj;(m£,^ 


)lk ■ 




(B35) 



Numerically we follow the following procedure : 

(a) Define the anomalous dimensions in Eqs. HAfflAfSp . 

(b) Define (71)} from Eqs. (|B22IB29p . 

(c) Define Eq. (|B30IB35^ . 

(d) Plug (a,b,c) into Eqs. p6IB21|) . 

This is much simpler than inserting the explicit formulae of Appendix lUl below. 



APPENDIX C: EXPLICIT RGES FOR THE SOFT SUPERSYMMETRIC BREAKING TERMS 

T he ex plicit RGEs for the soft supersymmetric breaking terms have appeared also before in Refs. 8] and |l07j . 
Ref. |l07j contains the full set (aside from the aforementioned S term), but we disagree with several terms in the 
equations for (m^ ^ ) and (mg^)ij. Ref. Q is restricted to contributions of the third generation quarks and leptons. 
We arrange here the explicit formulae of the full (not flavour dominance assumed) RGEs. As a cross check, we 
have calculated them by first using the explicit formulae of Ref. and second by using the method described 
in Appendix ^ We found agreement using both methods. Thus the RGE for the bilinear /i and k,; terms of the 
superpotential parameters is given by 



lOTT — = n 
dt ^ 



3 Tt{YuY\j) + Tr(3Yz5Yl3 + YbY^,) - -gl - 35: 



AJ,„Yb + 3A^„Yz, 



(CI) 



pn 



IGtt 



2 dt^i^ 



3 Tt{YuYI) ~ -gl - 3.9: 



Y£;Y^ + A£„A^„ +3Ac„a1 



A£„Y|; + 3A^„YB 



(C2) 



Similarly, the RGEs for the soft SUSY breaking bilinear terms can be read from. 



167r2 — = B 
dt 



3 Tr(Yt,Y£/) + 3 Tt{yIYd) + Tr(YljYB) - -gf - 3gl 



+ M 

- Di 



6 Tr(Yj,h[;) + 6 Tr{Y\jhn) + 2 Tr(Y^hB) + -gfMi + QglM2 



Yb + 3Al5„ Yc 



Hi 



In 



2A^„hE + 6Al,„hfl 



,dD. 



dt 



3Tr{YuYl)-^gf-3gl 



+ Hi 



6 Tr(hc;Yj,) + -gfMi + 6^1 M2 



- 2/x 



YbY^ + Ab„ A^„ + 3Ad^aI„ 



+ 2ki 



hfiY^ + hE„A^„ + 3hj3"A|5„ 



hEnY^+3hDnY^5 



A^n Y^ + 3A£)nY^ 



The RGEs for the soft SUSY triUnear couphngs are given by 
2d{hE)ij 



dt 



167r 



2d( h£))jj 



167r 



2 



+ 
+ 
+ 
+ 
+ 
+ 



+ 
+ 
+ 
+ 
+ 
+ 



J^Ei ) 

YeYI + Ae-AI^+3Ad-A^d 



hfij )ii 

YE)il 

^E)ij 
^Ei)il 

^D)lj 
hD)ij 



Tr{Y],YE) + 3 Tr(Y+,Yc) - -gf - 3g. 



In 



-A*e„Ye - 3A^„Yd 
4(Yt,hB)y+2Tr(At,, 
2hBY; + 2hBnA;„ + 6hD.At,„ 



2 Tr(Y^hB) + 6 Tr(Y]3hz5) + + 652^^2 



-2(A|;„hB) - 6{AlnhD) 



In 



2{Y\jYDh + 2 Tr{Al,An,) + 2{AunAlj^)ji 



YdYI+YuYI 



Al„ADn 



16 



Tr(Yt,YB) + 3Tr(Yt,Yo)-^52_3^2 ^ 



(A^„Yb)+3(AJ,„Y^) 



In 



Yn)uU{Y\,hD)ij + 4 Tr(Aj,,h^,) + 4(hc;„A^„)iJ 



Ani)ii 



2hz5Y^ + 2huYy 



2A^„hcn 



14 
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2 TrCYljhij) + 6 Tr(Yt,hz5) + j^ffi^i + 65IM2 + y^l^a 



2(A^„hB) + 6(AJ,„hz5 
2(Yt,Yt,)y + Tr(At,, 



167r- 



dt 



1677 



2^( hj)fc)ij 
dt 



+ 
+ 
+ 
+ 
+ 



+ 
+ 
+ 
+ 



+ 
+ 
+ 
+ 



+ 
+ 
+ 
+ 
+ 
+ 



hc/)y 
^u)ij 



3 Tr(Yt,Y^) - ^^gl - 3gi - ^-^gj 



26 S2 
6 TiiYljhu) + -5?Mi + 6glM2 + —gll^h 



KEnY% + SAd'^ Y^, 

- in 

YeYI - A£.A^„ - 3Az3nAt,„ 



YbY^ + KEnK^En + SA^n A^j^ 



— A^n Y^ — SA/JnY^ 

4(Yt,h£),fc+2Tr(A;,hB0 

2hEnY|; + eh^n Y^ 

-2\iEY^E - 2hBnA^„ - 6hz3.At,, 

2hEY^ + 2hBn AJj, + 6hz5n AJ,, 



-2hEnY^ - 6hcn Y^ 



r5f + 3^2 



+ 2Tr(A^,AB.) + 2(A^,„ 
YbY^ + Ajsn a]j„ + 3Ar,nAi,„ 

ABnY^ + 3ABnY^, 

(Y^yJ,),, + {YuY]j)^i + {A^j,^Ad^)i, 
4(Y|,ho);fc + 4 Tr(At„hB.) + 4(Aj,„ht;n) 
2h£;Yt, + 2h£. A^5„ + 6hi3.Aj, 

2hEn Y^ + 6hzjn YJ, 
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- {hjjk)ij 

+ (A£)fe)jj 



7 



dt 



+ (h[/')fei 



1 o 2 2 
Yg^l + 352 + -3-53 



2(Y|)Y£))/fc + 2 Tr(Aj,,A^.) + 2(A^„ 
2(Yt,Y[;)H+ Tr(At^,Ac,. 
-2(Yl,Yc)y - 2 Tr(Atj,A^. ) - 2(At,„ Ac/„)j, 
4(Y|,hB)ifc + 4 TV(A^,hj3.) + 4(Aj,^ha-) 
4(Yt,h[;)(, + 2Tr(A^,h[;.) 
-4(Y^ho);,- - 4 Tr(Aj,,hB.) - 4(A[,„han)i,- 



(C9) 



(Ac/')jfe 



(CIO) 



The RGEs for the gaugino masses are not affected by the ^p-coupUngs up to 1-loop. The RGEs for the SUSY soft 
breaking masses are given by, 

lg^2^(^]^Jii = 2(Yt,YE),„(mg2)„,- + Tr(At,,ABn)(mg')„,- 

+ 2(mg^)i„(Y^YB)„,- + (mg2)i„ TV(A^„Ab^) 

+ 4(Y^YE)i,ml,^ +4(Yt,A^,),,(m^^^^) 

+ 4Tr[(m£')At,,AB,]+4(m|^^J(At,, 

+ 4[Y^(m£2)Yi,]i,-+4(ht.h + 2 Tr(4, 

= (m£2)i„(YsYl,)„,- - (m|^^J(A^,Ys),-, 



(Cll) 



IStt 



2^( mL )i 
dt 



In] 



+ (mL^)i„(AB,Ajj,)„j + {Y EY'^E)in{l^t 

- (Ai5,Y^)i9(m^^^.) + (AB<,A^,)i„(mL^)„j 
+ 2(YsYt,),,.ml,^ +2(AB.)i.(Yt,)p,.(m^^£^) 

- 3(Az,,Y^,)i,(m^^£.) +3((m£')Az,,At,,),,. 

- 3(AJ,<,Y^),g(m|^^^) + 3(A,,,At,,(mL')),,- 

+ 2(Yi<;)ip(A^^) 

+ 2(Yi<;),.(mg'),,(Yt,),,. + 2(AB0ip(mE').5(At 

2\ 



ipV"'E IrqV^Ei IP] 
2\ 



+ 6(A£)r)jp(mg^)rg(A|j,)pj- +6(A£,fc)ji(mQ^)„;(A|jfe)„j 

+ 2(hBh|j)i,- + 2(hB<,h^<,)i,- + 6(hz3<,hJ„)y 



(C12) 
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+ 3{An.AUm{ml^^J + [ Tr(Y^,Yt,) + 3 Tr(Yt,Yz,)](m^^^^) 

+ (mL^)™(A^<, Yb - 3Al),YD)nq + 2(A^p)jq(mL^)gr(Y£;)rp 

+ 2(A^,Y£;)ir(mg^)rg - 6(A|),Y£,)ir(m£,^)rq 

- 6(A^p)i,(mQ^)gr(Yr,)rp - [2h|;,h£; + 6h^,hi5]i, , (C13) 

IGtt^ J ' = [Y,5Y|,+Yc;Yt,],y(mQ'),„ + (At,,Ai,,)jn(mQ')i„ 

+ [YdYd + Yc/Y^]i„(mq^)„j + (mQ^)„j (Ajj^Ar,,)^^ 

+ 2(Yz,),,(mD').,(Yl,),, +2(Y,3Yl3),,m|,^ 

+ 2{A^jjp)jq{mf^'^)qr{ADp)ri + 2{Yu) 

+ 2(Yc/Yj,)yTO2^^ +2(A^,Az3.)j.(mD')., 

- 2(m^^£^)(A^.),,(Yt,)fe, - 2(Yz,).,(At,,),fe(m|^^^) 

+ 2[hr,h|) + huhlj]^.j + 2(h^,hi5,)ii 

- (^5?|Mip + 6gl\M,\' + §glm' - IgfS^S.j , (CM) 
16^2 <mD^)ii ^ 2(Yt,Yc),„(mg^)„,+2TV(At,,Az,n)(mg2)„, 



+ 2(Ac/.A[„)„i(m5')„,- +2(yJ,Yb)„,( 

IHfi Jin 

+ 2 Tr(ABj A^„)(mfi2)i„ + 2{Aup 

+ <'y]D)i<li^Cl\r{YD)rj+^{Y\j)ip{YD)pjrn%^ 

+ 4(A^,Aj3,),r(mQ^) 

qr 

+ 4(A|^p)gi(Ac/p)j>(mg,^)^, + 4(A,yiA^,)ji(m^^)i, 
+ 4(hl,hi5)»j + 4 Tr(h^,hi),) + 4(hc/ph[,p),, 

- A{AyYnhinil^i^)-4{AD.Y*Mrnl^i^) 

- (^ff?|Mip + §9i\Ms\' - lafSy^ , (C15) 
^Q^2d(ni^ = 2(Yt,Y^)„,(mu')„, + Tr(At,. A^„)(mu^)„, 



2(Y^Yt;)„,(mu')„ + Tr(At,.A[,„)(mu').„ 
4(Y^,),,(mQ'),,(Yc;)., 

4(A[,.At;.),,(m6'),g + 4(hJ;h[;),, + 2 Tr(h^.hy.) 



|5?|Mir + |^?f|M3r + ^5?5)<5,,, (C16) 

^g^2^ = 2 Tr(Y^Yj,)mlf^ + 6 Tr(Yt,Yc)ml,^ + (Y^AfiO^^m^^iJ 
- 3(A^.Y^,),fe(m^^^_^) - 3(A^.Yz,)gfc(m|^^^) 
+ (A^Y ) + 2(Yt,Y 
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where 



and 



+ 6(Yz5Yl,),,,(mQ2),, + 2 Tr(h^h£;) + 6 Tr(ht,hc) 



^5?|Mip + 6g||M2p + |9?5) , (C17) 
= 6Tr(Y^Yt;)m|,^ +6(Y^,Yt,),,(mu'),, 



6(Yc/Y^,),,(mQ'),, + 6 Tr(ht,h[;) 



^5?|Mip + 6g||M2p-^g25) , (C18) 



5 = m^^ - m^^ + Tr[mQ2 _ rn^^ _ 2mu2 + rng,^ + mg^] . (C20) 



APPENDIX D: FOUR-BODY r-DECAY 



In this Appendix we compute the four-body decay — )■ r^/^+ud via the _^p-operator X' L^^QiDi. The relevant 
Feynman diagra ms a re given in Fig. |Sla,b,c. We neglect the contributions from the heavier neutralinos. Using the 
notation of Ref . |l08l | the three amplitudes corresponding to Fig. |S1 are given by: 

= + (^2_ jp|^'2_M2) ('^'^^") {^[^^^^ + + ^Wl^] (Dl) 

■^^ = -7T^ '^-TU^^^^7TA^PL^^){^{^rPL + brPR}{i + M■^)PLd\ (D3) 



Here the four-momenta are denoted by the particle symbol. The momenta jl, u, d, x flow along the corresponding 
propagators from left to right. We use below that x = /Lt-|-M-|-das well as the notation Np = — for the 
denominators in the propagators. We have assumed there is no mixing in the scalar /i, u, and d sectors. However, we 
allow for mixing in the stau sector, /i r, ttf,, dp a re the only sparticles that couple to the R-parity violating operators. 



fci r% 

The coupling constants are given by |3, ISo, Il08l | 



COSt^w / \ COS0W / 



bf, = -eTVii — , bu = ee„7Vii H ^ , (D5) 

COS t^w cos dw 

gcd sin^ 9^N{1 

ad = -eediVii H . (D6) 

cos t'w 

The total matrix element squared is given by 

\M\^ = [ \Ma\^ + \Ma\^ + \Ma\^ + 2 3? (MaMl + MaMl + MbMl^] , (D7) 
where A'c = 3 is the colour factor and 

-|C\/2|L 12 

\Ma\' = ^^^|^d.^[|a.|X^.^+ |6.|Vr,X,/^,x)] , (D8) 
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(c) 

FIG. 8: Feynman Diagrams for the decay f ^ TixiY ^ T{pud) via the operator L^QiDi. 



i(3\/2|L 12 

l-^^l' = ;v2Af2 + \br\'g{r,X,u,x)] , (D9) 

= [|a.pM^r.d + \br\'giT,x,d,x)] , (DIO) 

X d 



25R 



(M.Mi) = - jV^'^r [\ar\^M^g{T,fi,d,u) + \br\' fiT,x, ^^,d,u,x)] , (DU) 
25R(M,Mt) = ^ y [\ar\Hd'^g{T,fi,u,d) + |&.| X, M, x)] , (D12) 

X d 

23? (Al,Alt) = ^^^2]^ [kr pM^5(r, /X, d) + |6.| V(t, X, x)] • (D13) 

The functions are given by 

g{a,b,c,d) = a-b c-d - a-c b-d + a-d b-c , f{T,x,a,b,c,x)^-x'^g{T,a,b,c) + 2T-xg{x,a,b,c). (D14) 

The squared ampUtude in Eq. ljD7|) can be used in Monte Carlo simulation programs to generate events with a decaying 
stau. We are here interested in an analytic approximation for the total decay width. To this end, we shall assume 
X^ «C . This is equivalent above to setting 6^ = 0. Furthermore we assume that all scalar propagators are 
dominated by their mass terms and the scalar fermion mass is universal: to^ = rriu = = rh. In this simplified case 
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the amplitude squared is given by 

2 _ 16A'2|a^|2iVc 



\M\' 



-btj.bl g{T, /i, d, u) + hf,a*^ g{T, /i, u, d) + b^a*^ g{T, u, ^, d)] 



The total width is given by [109l lllO| 



(D15) 



(D16) 



where fci = t, A;2 = ^3 ^ u, = d. After the simplification our matrix element squared consists of three kinds of 
terms which depend on the final state four-momenta: (r ■ ^){u- d), (r • • d), and (r • • m). As can be seen from 
the phase space integral, these all contribute the sameLthey simply correspond to a relabeling. We thus explicitly 
integrate only the first term. Using Eq (4) from Ref. ,109] with iV = f — fci — /c2, we see that 



/ 



and we thus obtain 



n 



d^h 

2^ 



S^{N -u-d) [t- n) (u-d) 



d^T r d^ 



2Er J 2E 



In the rest-frame of the decaying stau with the z-axis in the direction of the 3-momentum of the r 

T = (Aff, 0,0,0), T = 0,0,1), ^i = i;^(l,sin6i,0,cos6') . 

Performing the integrals over d^r and d(j)^ 

Ai = Y J dErJ dEf,J dcos 9 E^E^il - cos 9) [M? - 2MfEr - 2MfE^ + 2E^Er{l - cos 9)] 



(D17) 



(D18) 



(D19) 



(D20) 



It is convenient to change to dimension less v ariables E'^ = ^Mf^;, Er = ^Mfy, and 1 — cos 61 = 2w |109| . Implementing 
the integral boundaries given in Refs. [Toflirnnj . this leads to the result 



Ai 



dz 



dy / dw + dy 



dw 



1-2 J {y+z-l)/yz 



[z'^y'^w{l — z — y + yzw)\ 



^3 Ml 



X .(D21) 

25 720 ^ ' 



We thus have for the total width 



257r5M2rfi4 



(D22) 



where K = 1/(720 x 2^) = 1/23040. 
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